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ABSTRACT 


A theory  la  developed  whl.ch  describes  the  behavior  of  a 
vapor  bubble  in  a liquid.  Its  physical  basis  la  the  assumption  that 
the  heat  transfer  effects  which  accompany  the  evaporation  occurring  at 
the  bubble  wall  when  the  bubble  grows » or  the  condensation  that  occurs 
there  when  the  bubble  collapses,  are  dynamically  important.  The 
basic  equations  of  hydrodynamics  are  shown  to  reduce,  for  the  problem 
under  consideration,  to  a dynamic  equation  which  describes  the  behavior 
of  the  bubble  wall,  and  a heat  convection  equation  for  the  liquid  which 
is  crupled  to  the  dynamic  equation  by  a boundary  condition  at  the  bubble 
surface.  A solution  for  the  heat  problem  is  obtained  under  the  assumption 
that  significant  temperature  variation  in  the  liquid  occurs  only  in  a 
thin  thermal  boundory  layer  surrounding  the  bubble  vail.  An  estimate 
of  the  correction  to  the  temperature  solution  is  also  derived.  Once 
the  temperature  at  the  bubble  wall  is  given,  the  vapor  pressure  within 
the  bubble  is  known  and  the  dynamic  problem  becomes  detsrminatci . 

The  theay  la  applied  to  the  eases  of  the  growth  of  a vapcr 
bubble  in  a superheated  liquid,  and  the  collapse  of  a vapor  bL..blo  in 
a liquid  below  its  boiling  temperature  at  the  exVjrnal  pressiu'e.  The 
sliqpllX^^ng  physical  assumptions  made  in  the  course  of  the  Invoatigation 
ai^  Justified  for  the  specific  example  of  vapcar  bubble  behavioi*  in  water. 

A comparison  of  the  theory  with  experiment  is  given  for  the 
observable  range  of  bubble  growth  in  superheated  water,  end  the 
agreement  is  found  to  be  very  good. 
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I.  INTRODUCTION 


The  term  cavitation  is  used  to  describe  the  presence  of  a vapor 
phase  in  a region  filled  predominantly  with  liquid.  The  condition  necessary 
for  the  appearance  of  cavitation  is  that  locally  the  vapor  pressure  of  the 
liquid  must  exceed  the  external  pressure  to  which  the  liquid  is  subjected; 
this  condition  is  by  no  means  sufficient  to  produce  cavitation,  however, 
because  of  the  presence  of  s\arface  stresses  in  the  liquid.  These  stresses, 
attributable  to  the  short  range  attraction  of  the  liquid  molecules  for  one 
another,  tend  to  reduce  any  surface  element  of  the  liquid  to  one  having  the 
least  (mean)  curvature  consistent  with  the  mechanical  constraints  imposed 
on  the  liquid.  Thus,  an  otherwise  unconstrained  vapor  cavity  will  be  spherical* 
The  resultant  of  the  stresses  on  an  element  of  surface  is  a force  directed 
along  the  normal  drawn  from  the  concave  side  of  the  element.  For  a vapor 
cavity  to  grow,  the  vapor  pressure  must  compensate  not  only  the  external 
pressure  on  the  liquid,  but  also  the  effective  pressure  of  the  surface  stresses. 

Since  the  surface  stresses  increase  with  the  curvature  of  the 
surface,  there  is  a minimm  possible  size  for  an  imconvStrained  pure  vapor 
bubble  existing  in  the  liquid,  even  at  temperatures  above  the  boiling  point 
of  the  liquid  at  the  prevailing  external  pressure.  Smaller  bubbles  are 
unstable  against  collapse.  The  question  therefore  arises  as  to  how  a bubble 
could  form  initially.  The  problem  of  the  nucleution  of  vapor  bubbles  has 
been  extensively  studied  in  recent  years,  notably  by  Harvey and  Pease, 
in  connection  with  their  research  in  animal  physiology.  The  conclusion 
drawn  from  those  studies  is  that  in  a moderately  superheated  liquid,  the 
nuclei  for  cavitation  bubbles  consist  of  small  permanent  gas  bubbles  in  the 
liquid,  or  gas  pockets  stabilized  on  solid  purticloe.  When  those  are 
removed  from  the  liquid  (by  agitation,  continued  boiling,  or  by  compressjog 
the  liquid  under  pressures  grout  enough  to  forco  the  gases  into  solution) 
cavitation  ceases,  and  cun  bo  reintroduced  only  by  subjecting  the  liquid 
to  extreme  tension  or  high  temperature.  Thus,  water  put  under  u pressure 


of  several  hundred  atmospheres  for  a period  of  a few  hours  boeoraeu  able  to 

(3) 

withstand  negative  prossuroo  an  great  as  150  atjn.  without  rupturing,^  or 


o ) 

can  be  heated  to  270  1 before  it  explodes,' 


The  residual  nuclei,  following 


degassing,  are  believed  to  consist  of  hy'lrophoblc  aubstuncoo  in  the  liquid 

(2) 

or  at  Ita  surfaoo.' 
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The  process  involved  in  the  nucleation  of  a bubble  by  a permanent 
gas  can  be  explained  by  a simple  model.  Suppose  the  gas  satisfies  the 
perfect  gas  law.  Then  the  pressure  p of  the  gas  in  the  bubble  is  given 
in  terms  of  the  temperature  T and  radius  R of  the  bubble  by 


(1) 


where  N is  a constant,  proportional  to  the  n\amber  of  moles  of  gas  in  the 
bubble.  The  pressure  p^  due  to  surface  tension  is 


P 


s 


(2) 


being  the  surface  tension  constant  of  the  liq\iid.  If  the  bubble  is  in 

equilibrium,  the  vapor  pressure  is  a function  of  T alone,  say  p (T). 

eq 

Denoting  the  external  pressure  by  p^  , one  has  as  the  condition  of 
equllibrivim  that 


h ^ = Pa  Pc.  • (3’ 

The  equilibrium  will  be  stable  if  the  pressure  difference  Pg+  Pgq"  Pg“  Poo 
is  a decreasing  function  of  the  bubble  radius  at  the  point  of  equilibrium. 
These  conditions  are  conveniently  expressed  in  terms  of  a function 

fjCR)  = ^ [ P«-^  Ps  - PsqCr)}  = T [ Cp„-  Peq(^)l  + Zo'/j.  (4) 

Thus,  for  a given  gas  content  N,  the  equilibrium,  radius  (or  radii J 
of  the  bubble  is  given  by 

f^(R^)  = N,  (5) 

and  the  condition  for  stability  becomes 

^ Of 

wh0.ro  f|(R)  denotes  the  derivative  of  f,p(R)  with  respect  to  R. 

Bolow  the  boiling  point  of  the  liquid,  - p,  > 0,  and  so  f,j,{R) 
is  an  IncroQOing  function  of  R.  Thus,  there  is  just  one  equilibrium  raditm 


= 3 “ 


of  a gas  nucleated  bubble  for  a given  value  of  T below  the  boiling  point 
of  the  liquid t according  to  {5)»  and  the  equilibrium  is  stable  at  that 
radius  by  (6)«  Above  the  boiling  point  of  the  llquidt  the  coefficient  of 
in  fj(R)  is  negative.  Hence  as  R increases  for  fixed  T,  f^{R) 
increases » reaches  a maximum^  then  decreases.  Accordingly i eq.  (5)  may 
give  two  equilibria  radii  > the  larger  corresponding  to  unstable  equilibrium » 
or  oie  equilibrium  radius  vdiioh  is  stable'  against  collapse  but  unstable  for 
gfbwthi,  or  it  may  afford  no  equilibrium  radius.  Inasmuch  as  Pgq(T)  is 
an  increasing  function  of  T,  f^(R)  is  a decreasing  function  of  T for 
ai^  fixed  R|  so  that  the  curves  of  f^(R)  on  an  f - R diagram  form  a 
[hotiintersec  ting  family  except  for  the  common  point  at  the  origin  * and  the 
'Curves  for  large  T fall  below  those  for  small  T«  In  particular*  the 
liimiximum  of  f^(R)*  which  occurs  when  the  liquid  is  heated  above  its  boiling 
point  at  the  external  pressure  p^  » decreases  with  an  increase  of  T.  A 
typical  f - R diagram,  drawn  for  water  at  1 atm.  external  pressure,  has  been 
presented  in  Tig.  1 to  illustrate  these  general  remarks i 

Consider  a gas  nucleated  bubble  which  is  in  stable  equilibrium 
ip  a liquid  below  its  boiling  point,  and  suppose  the  temperature  to  rise 
slowly.  The  bubble  radius  will  then  increase  steadily » with  the  bubble 
remaining  in  stable  equilibrium  as  the  temperature  increases  past  the  boiling 
point*  until  there  is  finally . reached  a critical  temperature,  and  a 
corresponding  critical  radius,  above  which  the  bubble  cannot  es^st  In  stable 
equilibrium.  A further  increase  in  tenqjerature  releases  the  bubble  for 
dynamic  growth. 

ThS  nucleation  process  described  can  be  understood  on  the  basis 
of  the  f - R diagram  of  Tig.  1.  The  locus  of  the  process  for  any  given 
bubble  is  a horizontal  line,  whose  ordinate  is  fixed  by  the  gas  content 
of  the  bubble.  At  the  beginning  of  the  process,  the  bubble  is  represented 
by  the  intersection  of  the  given  horizontal  line  with  the  f^CR)  curve 
for  the  initial  temperature.  As  the  temperature  increases,  the  point 
representing  the  lAibble  shifts  to  the  right  on  the  f - H diagroa  to  f^(R) 
curves  drawn  for  higher  temperatures,  and  correspondingly  the  bubble  radius 
increases.  The  process  terminates  whan  the  bubble  point  reaches  the  f|j.(R) 
curve  has  o ndxinum  at  the  ordinate  of  the  horlzcmtal  lino.  (The  locus 
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- f-R  diagram  for  water  at  I atm.  external  pressure. 
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of  naxima  is  represented  In  tha  diagram  by  a dashed  lias , the  "locus  of 
critical  points".)  The  bubble  considered  is  then  stable  against  collapse 
at  the  temperature  of  this  critical  f^(R)  curve  and  at  the  radius 
of  its  maximum  point,  but  is  in  a situaticai  of  unstable  equilibrium  vdth 
respect  to  growth.  A further  increase  in  ten^jerature  upsets  the  equilibrium 
and  releases  the  bubble  for  dynamic  growth.  Since  the  surface  tension 
effects  relax  with  an  increase  of  bubble  radius,  the  bubble  (which  is  now  in 
superheated  liquid)  will  continue  to  grow  indefinitely. 

The  nucleation  process  is  qualitatively  similar,  but  from  an 
analytic  standpoint  more  conplex,  when  the  cavitation  bubble  grows  from 
a solid  particle  in  the  liqiiid,  or  from  a gas  pocket  stabilized  on  a solid 
particle.  The  bubble  may  pass  through  several  intermediate  stable  or  un- 
stable equilibria,  depending  on  the  size  and  shape  of  the  particle.  Since  a 
pure  gas  bubble  will  eventually  rise  because  of  gravitational  effects  and  so 
be  removed  from  the  liquid  as  a source  of  nucleation,  the  majority  of  cavita- 
tion bubbles  may  be  supposed  to  grow  from  nuclei  containing  solid  particles. 

The  mechanism  discussed  above  for  the  release  of  a cavitation  bubble 
for  dynamic  growth  is  the  counterpart  of  boiling.  One  may  analyze  in  a 
similar  manner  the  shift  of  a bubble  from  stable  to  unstable  equilibrium, 
and  its  release  for  growth,  by  a decrease  in  the  external  pressure.  An 
equivalent  process  occurs  in  cavitating  liquid  flow,  the  pressure  drop  in  the 
vicinity  of  the  bubble  nucleus  being  caused  by  a change  in  the  flow  pattern 
due  to  the  presence  of  a submerged  obstacle.  In  this  case,  however,  the 
bubble  does  not  ordinarily  continue  to  grow,  but  is  forced  to  collapse  by  a 
pressure  rise  which  follows  along  the  path  of  the  bubble.* 

A different  phase  of  the  nucleation  problem  has  been  investigated  by 
Glaser, who  used  degassed  diethyl  ether  at  1 atm,,  superheated  100°C 
above  its  normal  boiling  point  of  34°C,  as  the  working  fluid  in  a Bubble 

Chamber  designed  to  locate  the  path  of  a charged  atomic  particle.  A series  of 

^ - - __  . ..  . 

Considerable  local  pressures  can  develop  at  the  point  of  bubble  collapse. 
If  the  bubble  collapses  near  the  surface  of  a submerged  object,  the  sudden 
vuibalunce  of  pressure  resulting  may  be  sufilcient  to  dislodge  particles  from 
the  surface.  For  a recent  study  of  cavitation  damage,  see  M.S.  Plesset 
and  A,  Ellis,  Proceedings,  Annual  Meeting  ASME,  December  1954. 
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vapor  bubbles  appears  along  the  traolf  of  gueh  a particle  in  the  liquid*  If 
photographed  a few  nioroseoonds  after  the  detection  of  the  particle » the 
track  is  fairly  well  defined  by  the  babbles.  The  physical  aechanlsm  of  the 
bubble  nucleation  in  the  Bubble  Chamber  has  not*  as  yet»  been  fully  explained. 

In  treating  the  problem  of  bubble  growth,  it  becomes  necessary  to 
make  some  assumption  concerning  the  nucleation  process.  The  simplest 
assmptlon  to  moke  about  the  bubble,  and  the  one  which  will  be  made  in  the 
analysis  to  follow,  is  that  the  bubble  contains  initially  no  permanent  gas 
or  Solid  particle  nucleus.  It,  of  course,  follows  from  this  assumption 
(N  = o)  and  eqs.  (4)  and  (5)  above  that  there  is  no  radius  of  stable 
equilibrium  for  such  a bubble*  Nevertheless,  if  the  liquid  is  heated  above 
its  boiling  point  there  will  still  be  a radius  of  unstable  equilibrium, 
which  satisfies  the  eqixation 

f;j.(Ro)  = 0.  (7) 

The  radius  of  the  pure  vapor  biibble  given  by  eq.  (7)  is  related  to  the 
critical  radius  R^  for  Tinstable  equilibrium  of  a gas  filled  bubble  at  the 
' Same  superheat  by 


While  it  is  not  possible,  physically,  to  form  a pure  vapor  bubble  at  the 
radius  R^,  the  details  of  growth  of  such  a conceptual  bubble  differ  in  no 
essential  way  from  those  of  a gas  filled  bubble  growing  from  unstable 
equilibrium,  or  from  those  of  a bubble  otheridse  nucleated.  Equilibrium 
bubble  radii  for  a pure  vapor  bubble  in  water  at  1 atm.  external  pressure, 
as  a ftinction  of  the  water  temperature,  will  be  fovmd  in  Table  II  (p./65) 

The  equllibri-um  radius  of  the  pure  vapor  bubble  and  the  corresponding  oritical 
radius  of  the  gas  nucleated  bubble  are  decreasing  functions  of  temperature 
(see  Hg.  1). 

The  method  considered  in  the  analysis  below  for  the  release  of  an 
equilibrium  vapor  bubble  for  growth  will  be  an  increase  cf  liquid  temperature, 
rather  than  a decrease  of  pressure.  One  reason  for  this  choice  is  that  the 
tenqwrature  change  involved  in  the  heating  of  a liquid  usually  proceeds  at  a 
sufficiently  low  rate  that  it  ceases  to  influence  the  behavior  of  a bubble 
sqon  after  the  growth  of  the  bubble  has  begun,  ihe  bubble  growth  then  becomes 


§9*L?^eterH!,lnsd  j **#!!?ss*s  When  the  b'ibble  growth  Is  Isltietod  ^ 

parasaure  ebangest  this  may  not  be  the  ease. 

Another  reason  for  the  oh<^ee  is  the  availability  of  experimental 
data  for  the  growth  of  cavitation  babbles  in  superheated  water « with  which 
the  predlotlona  of  the  theory  developed  below  may  be  ampsred  (see  ilg8«  4t 
5$  6)«  After  the  growth  of  a oavltatlon  bubble  has  beguni  the  details  of 
nueleatlon  become  unln^ortant*  The  bubble  tends  to  become  spherical » and 
is  adecpmtely  represented  by  the  pure  ^por.  bubble  model  used  here* 

The  process  of  growth  of  a ^vitatlcn  bubble  in  a superheated 
liquid  may  be  described  as  follows i Uhen  tlie  bubble  is  at  its  critical 
radius  (R.  for  a gas  nucleated  bubble*  the  radius  R_  of  eq»  (7)  for  a 
pure  vapor  bubble)*  it  la  tmstable  against  esgpanslon*  and  a sli^t  tea^esra- 
ture  increase  will  alert  the  bubble  growing*  The  initial  phase  of  growth 
is  oharaoteriaed  \if  the  relaxation  of  the  effective  pressure  due  to  surface 
tension  with  an  increase  of  bubble  radius*  The  pressure  unbalance  causing 
the  bubble  growth  is  thereby  increased*  end  correspondingly  the  rate  of 
expansion  increases  rapidly*  In  order  for  the  bubble  to  grow*  however* 
evaporation  must  take  plaoe  at  the  bubble  wall*  Because  of  the  latent  heat 
requirement  of  evaporation,  this  requires  the  temperature  at  the  bubble 
Tfall  to  drop  below  that  of  the  bulk  liquid*  which  in  turn  decreases  the 
vapor  pressure  at  the  bubble  surface*  Whether  or  not  the  decrease  in 
pressure  causes  the  velocity  of  the  bubble  wall  eventually  to  decrease 
depends  upon  the  rate  of  Inorease  of  bubble  surface  area*  It  will  be  eham 
that  such  an  effect  occurs*  The  babble  radius  ultimately  becomes  propor-* 
tiooal  to  the  square  root  of  the  time  of  growth*  In  this  asyoptotio  range 
of  bubble  expansion*  the  tenperatore  at  the  bubble  wall  approaches  the 
boiling  point  of  the  liquid  at  the  ext^iv^al  pressure,  and  the  pressure 
difference  producing  tiie  bubble  growth  tends  to  zero  with  the  radial 
velocity  of  the  bubble  wall* 

If  the  bubble  growth  is  arrested  and  the  bubble  forced  to  collapse 
by  a sudden  inorease  in  the  exterior  pressure*  the  flow  of  vapor  and  the 
flow  of  heat  at  the  bubble  wall  are  reversed*  tkmdsnsation  of  vapor  at  the 
surtlaoe  of  ths  bubble  raises  the  temperature  thex^*  resulting  in  an  increase 
of  vapor  pressure  which  tends  to  slow  down  the  rate  of  coHapse* 


It  i5  thus  apparent  that  coupled  vith  the  dynamic  problem  there 
is  a problem  of  heat  transfer  between  the  liquid  and  vapor  \jhich  arises  when 
the  bubble  changes  size.  The  heat  problem  \d.ll  be  solved  approximately 
under  the  assumption  that  significant  heat  transfer  occurs  in  the  liquid  only 
in  a thin  shell  surrounding  th©  bubble  wall.  The  solution  is  presented  in 
section  III,  along  with  an  estimate  for  the  first  order  correction.  The 
assumption  of  a thin  thermial  boundary  layer  in  the  liquid  is  reasonable  if  the 
thermal  diffasivlty  of  the  liquid  is  sufficiently  low. 

Insofar  as  the  liquid  is  concerned,  the  bubble  grows  or  collapses 

because  of  pressure  variations  at  the  bubble  wall,  and  possibly  at  the 

external  surface  of  the  liquid,  which  set  the  liquid  in  motion.  Thus  the 

heat  transfer  problem  Involves  convection  effects,  A treatment  of  the  heat 

transfer  problem  which  neglects  convection  has  been  given  by  Forster  and 
(7) 

Zuber,'  who  use  the  model  of  a stationary  liquid  containing  a moving  heat 
source  (coi»responding  to  the  moving  bubble  wall).  The  diffusion  solution 
for  the  heat  problem  obtained  from  this  model  leads  to  unrealistic  predictions 
for  the  rate  of  bubble  growth.  An  analysis  of  the  diffusion  solution  is 
also  presented  in  section  III, 

The  dynamic  problems  considered  here  are  the  growth  from  wstable 
equilibrium  of  a pure  vapor  bubble  in  a superheated  liquid,  and  the  collapse 
of  a vhpor  bubble  in  a liquid  below  its  boiling  point.  The  bubble  which 
collapses  Is  taken  to  be  at  rest  initially,  and  in  this  respect  is  a model 
for  a cavitation  bubble  whose  growth  has  been  arrested  by  an  increase  in  the 
external  pressure  on  the  liquid.  The  model  used  differs  from  an  actual 
cavitation  bubble  in  that  the  liquid  temperature  is  assTomed  to  be  'oniform 
Vdien  the  collapse  starts.  The  temperature  field  in  the  liquid  for  an  actual 
bubble  depends  on  the  past  history  of  the  bubble,  and  if  non-uniform  will 
affect  the  Initial  period  of  bubble  collapse.  The  solutions  for  the  growth 
and  collapse  of  vapor  bubbles  are  presented  in  section  IV,  together  with 
experimental  verification  for  the  case  of  bubble  growth  in  supei'heated  water. 

For  the  quantitative  solution  of  the  heat  problem  and  the  dynamic 
problem,  several  simplifying  physical  assumptions  have  been  made.  The 
arguments  for  the  validity  of  the  general  assumptions  are  justified,  as  they 
appear  in  the  text,  for  the  case  of  cavitation  bubbles  in  water.  Two  basic 
assumptions  may  be  mentioned  here,  however^  these  are  that  the  motion  possesses 


- 9 “ 


spherical  symmetry,  and  that  the  motion  is  irrotational.  The  latter 
assumption  is  independent  of  the  former  since,  for  example,  it  is  possible 
for  a swirling,  eddy  type  of  motion  to  occur  within  the  vapor  of  a spherical 
cavitation  bubble.  The  experimental  avldence  indicates  that  such  motion,  if 
it  occurs,  does  not  influence  the  bubble  behavior.  The  asstamption  of  spherical 
symmetry  is  more  serious.  This  requires  in  principle  that  the  asymmetric 
effect  of  gravity  upon  the  bubble  behavior  be  Ignored.  Actually,  the  rise 
of  a vapor  bubble  against  gravity  is  extremely  slow,  so  long  as  the  bubble 
is  small.  Thus,  for  water  superheated  by  about  2*^0,  no  great  error  is 
introduced  by  the  buoyant  force  provided  the  bubble  growth  is  not  followed 

-I 

beyond  a radius  of  the  order  of  10  cm,  which  is  much  greater  than  the 
equilibrium  radius  of  about  1 ,5  x lo"^  cm  for  the  102°  vapor  bubble  in  water. 
Bubbles  released  at  higher  superheats  grow  appreciably  faster  than  the 
102°  bubble,  and  so  are  relatively  much  larger  before  gravitational  effects 
become  important.  The  collapsing  bubble  has,  effectively,  no  time  to  rise 
against  gravity  before  its  collapse  is  completed. 

The  emphasis  in  the  following  treatment  is  laid  upon  the  physical, 
rather  than  the  mathematical  aspect  of  the  problem.  Thus,  a complete  table 
of  the  integrals  appearing  in  the  text  has  not  been  given,  although  a few  of 
the  more  obscure  integrals  are  evalaated  in  the  Appendix.  From  a mathematical 
standpoint,  however,  it  is  felt  that  the equation  for  the  growth  of  a vapor 
bubble  in  a superheated  liquid  may  be  of  some  interest,  inasmuch  as  it  offers 
a tractable  example  of  a nonlinear,  integro-diffsrential  equation. 
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II.  TORMULATION  OF  THE  PROBLEM 


Basic  Equations. 

In  terms  of  the  fluid  density » > the  (vector)  fluid  velocity  v, 

the  temperature  T,  the  pressure  p,  the  internal  energy  e per  unit  mass 
of  flviid»  the  stress  tensor  P»  the  thermal  conductivity  k,  the  coefficient 
of  viscosity  ^ f the  time  t,  and  the  heat  q generated  per  second,  per 
unit  volume  in  the  fluid  by  absorbed  radiation,  the  fundamental  Eulerian 
equations  which  describe  the  behavior  of  a fluid  (liquid  or  vapor)  are: 

The  equation  of  continuity 

If  + y . y)!,  = o.  (1) 


The  equation  of  motion  (with  the  neglect  of  external  body  forces, 
such  as  gravity) 

dv 


^ it  = 


The  heat  equation 


= P :W+7*kVT  + q. 


(2) 


(3) 


The  thermal  and  caloric  equations  of  state 


jO(P»  T),  8 = o(p,  T).  (4) 

For  a Newtonian  fluid,  the  stress  tensor  is  given  by* 

P = -pi  + [W  + (^)»  - I I(V  • v)l.  (5) 


The  notation  used  hero  is  essentially  that  of  Gibbs,  with  7 
denoting  the  gradient  operator.  The  symbol  V • denotes  the  divergence,  V x 
the  curl,  and  in  the  case  of  the  rate  of  strain  tensor  ^ is  the  vector 
gradient  (a  dyadic).  The  term  P : yv  in  (3)  represents  the  trace  of  the 
product  of  the  stress  and  rate  of  strain  tensors.  The  term  (’*JV)*  in 
eq.  (5)  is  the  transpose  of 

For  the  definition  of  the  stress  tensor  and  a derivation  of  eqs.  (1), 
(2),  (3),  see  rHlne-Thomson,  Theoretical  fivdrodynamics  (ifecmillan  and 
Co.,  Ltd.,  London,  1949)* 
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In  these  equations,  ^ denotes  the  total  derivative  vdth  respect  to  tin;®,  os 
coi^uted  in  a reference  frame  at  rest  in  the  fluid  element  under  consider- 
atlonj*  thus. 


J. 

dt  at 


+ V 


'7. 


(6) 


Physically,  eq.  (1)  expresses  the  conservation  of  mass.  Eq#  (2) 
relates  the  inertial  reaction  of  the  g3.ven  elementary  fluid  mass  to  the 
surface  stresses  acting  upon  it  (in  the  absence  of  external  forces).  Eq.  (3) 
is  essentially  the  statement  of  the  first  law  of  thermodynamics,  relating  the 
increase  of  internal  energy  of  the  mass  element  to  the  work  done  on  it  by 
its  surroundings  in  changing  Its  shape  and  size,  the  heat  energy  conducted 
into  it  across  Its  surface,  and  the  heat  generated  in  it  by  absorbed 
radiation. 


By  ths  use  of  standard  vector  and  tensor  identities,  the  terms  in 
(2)  and  (3)  involving  the  stress  tensor  defined  by  eq.  (5)  can  be  reduced 
to  the  forms 


V • P = -S?p  + [|  V(V  • x)  (V  X 2)1 » 

PsW  = -p(V  * v)  + (V  / v)^  + ^ (V  » v)^ 

+ 2v  • Cj  - Y X(v  ;<  :2)  - y(v  • y)1  ^ , 


) (7) 


provided  the  coefficient  of  viscosity  Is  considered  to  be  constant. 
Because  of  the  smallness  of  the  coefficient  of  viscosity  of  water  and  its 
vapor,  viscosity  effects  are  not  expected  to  play  an  important  role  in  the 
behavior  of  water  vapor  bubbles,  and  will  ultimately  be  ignored.  The 
viscosity  terms  have  been  retained,  however,  so  that  order  of  magnitude 
estimates  of  their  importance  may  be  made. 

When  the  assumption  is  made  that  the  fluid  flow  is  irrotational 


V X Y = 0 


(B) 


Other  symbols  for  the  total  (substantial,  convective,  particle) 
derivative  include  D/Dt,  and  a dot  placed  above  the  differentiated 
terms  q • 
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the  forms  (7)  further  reduce  to 

V • P s + I ^ 7(V  • v) , 

Pj  '^  = -p(7  > v)  + ~ I (V  * v)^  " 22  • 7(7  • v)  ^ 


> (9) 


If  the  further  assun^tion  is  made*  in  the  case  of  the  liquid,  of 
incompressibility  (y?  = constant),  eq*  (1)  gives 


and  eqs,  (9)  become  simply 


V « V = 0, 


P - -VPs 


P ! 7v  = ^7^  v^. 


(10) 


(11) 


It  will  be  observed,  on  substituting  eq»  (11)  into  eq.  (2),  that  the  viscosity 
terns  disappear  completely  from  the  eq\aatlon  of  motion  for  the  case  of 
irrotational  motion  of  an  incompressible  fluid.  If  the  fluid  is  viscous, 
the  motion  may  still  be  influenced  by  viscous  heat  generation,  howe-ver, 
provided  significant  heat  transfer  effects  take  place. 


The  Problem  In  the  Liquid. 

It  follows  from  the  assumption  that  the  liquid  motion  is  irrotational 
(eq.  (8))  that  there  exists  a velocity  potential  ^ throughout  the  liquid, 
such  that 

V = -70.  (12) 

Since,  further,  the  liquid  is  taken  to  be  incompressible  (eq.  (10)),  the 
velocity  potential  is  a solution  of  Laplace’s  equation 

0 = 0,  (13) 

The  spherically  symmetric  soliitlon  of  eq.  (13)  is  of  the  form 

0 = + B(t), 


(14) 
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\Ajere  r is  the  radial  coordinate  from  the  center  of  the  bubble.  The 
liquid  velocity  corresponding  to  eq.  (14)  is  purely  radial,  and  given  by 
eq.  (12)  as 

V = ^ . (15) 

r 

If  the  liquid  velocity  at  the  bubble  wall  is  denoted  by  v(R)>^  eq.  (15) 
gives  the  relation 

A(t)  = R^(t)  v(R),  (16) 

where  R(t)  is  the  radius  of  the  bubble  surface,  Eq,  (15)  thus  becomes 

2 

v(r,t)  - v(R),  (17) 

r 

and  if  the  velocity  potential  is  norinalized  to  zero  at  r = «»  , eq.  ('14 ) 
becomes 

= s!(U.yta?. . (18) 

The  equation  of  motion,  according  to  eqs,  (2),  (11),  is 

5v 

/?[-“  + V • yv]  = -yp,  (19) 

^ at 

The  vector  identity 

v>  (y  X v)  = “ 2 • 7V 

and  the  assmption  y X v = 0 give 

V * yv  = ^ yv^j 

hence,  by  eq.  (12),  the  equation  of  motion  (19)  may  be  written 

+ =-^p, 

from  which  a Bernoulli  relation 

" it  2 

follows  by  integration,  the  density  having  been  taken  constant.  From  eqs,  (17), 

v(R)  is,  in  general,  not  exactly  equal  to  the  radial  velocity  R 
of  the  bubble  wall  (because  of  evaporation  or  condensation  occurring  there.) 


(18)*  the  left  side  of  eq,  (20)  vanishos  at 
eq.  (20)  la  simply  the  external  preasure 


_ ^ t . / a ^ - 

r = i«  • ao  tnuTV  xii 

divided  by  the  liquid  density* 


£!q.  (20}  thus  becomes 


G(t)  =jj^  • 

M 1 

at  " 2 ^ “ /o 


(21) 


Because  the  density  of  the  liquid  is  assumed  constant » the  Internal 
energy  e can  be  a function  only  of  the  temperature  T.  Over  the  limited 
range  of  temperatures  which  will  be  considered  here,  the  internal  energy  may 
be  considered  tu  vary  linearly  with  the  temperature.  On  neglecting  any 
constant  internal  energies,  the  caloric  equation  of  state  therefore  reduces  to 


e T,  (22) 

where  c^  is  the  specific  heat  (at  constant  volume)  of  the  liquid.  The 
thermal  conductivity  k,  like  the  specific  heat,  will  be  taken  constant 
over  the  temperature  ranges  considered.  Eqs.  (3),  (11),  und  (22)  thus  combine 
to  give  for  the  internal  energy  equation  of  the  liquid 


[|r  + V • 7T1  - k^~T  + + q-  (23) 

y V tJ\j  ^ ^ 

2 2 

In  the  analysis  to  follow,  the  viscosity  term  ^ in  eq.  (23) 
will  be  neglected.  By  using  the  solutions  thus  obtained,  it  is  possible  'to 
estimate  the  contribution  of  this  term  to  the  total  rate  of  heat  generation 
per  unit  volume  of  liquid,  'Ihe  specification  (17)  for  v gives 


and  accordingly,  the  viscous  heat  generation  is  a ranximiua  at  the  bubble  VKiU, 
whei’e  it  amounts  . to 


per  second,  per  unit  volume.  Here  v(R)  represents  the  velocity  of  the 

liquid  at  the  bubble  wall,  which  may  be  approximated  by  the  radial  velocity 

• 

R of  the  bubble  wall  Itself.*  The  coefficient  of  viscosity  of  water  near 
* 

This  approadJaation,  idiich  is  plausible  physically,  will  be  shown  to  be 
accurate  to  about  1 part  in  1000, 
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bubble  growth  to  be  presented  below  gives  for  a bubble  growing  in  water  at 

1 atm.,  superheated  to  103°C,a  maxinaiia  radial  velocity  §,  = 32  cm/sec 

when  the  bubble  radius  is  about  R = 3 x lo"’’  cm  (see  Fig.  7).  Combining 

these  figures,  we  find  for  the  rate  of  viscous  heat  generation  at  this  time 

time  ,2 

R** 

~ 12?^  = 4 X 10^  erg/oe*sec  ~ 10“*^  cal/oc*sec. 


But  the  total  tengsorature  drop  at  the  bubble  wall  near  the  time  of  maximum 
radial  velocity  (see  RLg.  8)  is  about  10^  °c/sec,  corresponding  to  a heat 
loss  from  the  liquid  at  the  bubble  wall  at  this  time  of 


dt 


~ 10^  cal/cc*sec, 


dtie  essentially  to  the  evaporation  occurring  there.  The  viscous  heat  gener- 
ation drops  off  sharply  away  from  the  velocity  maximum.  Clearly,  viscosity 
plays  a negligible  part  in  determining  the  growth  of  vapor  bubbles  in  water, 
Eq,  (23)  will  therefore  be  written  for  solution  as 

/°v®  + V • 7Tl  = k7^T  4 q . (24) 


The  foregoing  development  of  the  equations  for  the  liquid  has 
been  based  essentially  on  the  assumption  that  the  liquid  is  incompressible. 
The  validity  of  this  assumption  depends  upon  the  ratio  of  the  velocities 
attained  by  the  liqviid  to  the  velocity  of  sound  in  the  liquid.  For  the 
Rowing  vapor  bubble,  the  maximum  velocity  at  the  bubble  wall  is  never 
larger  than  a few  meters  per  second  for  the  highest  superheats  considered, 
so  that  con?)ressibility  effects  in  the  liquid  may  be  safely  ignored, 

For  the  case  of  the  collapsing  bubble,  several  of  the  asstimptions 
made  above  may  fail  near  the  point  of  collapse,  if  the  solutions  are  carried 
that  far.  Thus,  the  temperature  at  the  bubble  wall,  which  is  initially  balow 
the  boiling  point,  rises  sharply  near  the  end  of  collapse,  possibly 
* '' 

A tabulation  of  physical  constants  will  be  found  in  Appendix  A. 
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approaching  the  critical  temperatwe  of  the  liquid.  The  parameters yO»  k, 

c^j  etc.  of  the  liquid  cannot  be  taken  constant,  of  coxirso,  over  such  an 

extreme  temperature  variation,  nca*  is  it  valid  to  consider  the  liquid 

incoii5>ressible . Perhaps  the  most  significant  error  in  the  ba.sic  assunqstions , 

for  the  case  of  the  collapsing  bubble,  lies  in  the  fact  that  the  spherical 

(8) 

shape  is  inherently  unstable  near  the  point  of  collapse . ' The  collapsing 
bubble  tends  (theoretically  and  ex^jerimentally)  to  shatter  before  collapse# 
For  these  reasons,  although  the  assxunptiona  made  above  will  be 
retained,  the  calculations  for  the  collapsing  bubble  will  be  carried  only 
far  enough  to  indicate  the  trend  of  behavior  of  the  physical  quantities 
involved. 

The  Problem  in  the  Vapor. 

In  the  case  of  the  vapor,  the  main  simplifying  assumptions  are 
related  to  the  smallness  of  the  vapor  density  in  comparison  with  that  of 
the  liquid.  Thus,  the  physical  effects  of  the  vapor  inertia  may  be  expected 
to  have  a negligible  influence  on  the  rate  of  bubble  .growth  or  collapse. 

It  will  be  Shown  that  the  vapor  may  safely  be  considered  to  be  in  a state 
of  thermal  and  dynamic  equilibria,  insofar  as  its  internal  behavior  is 

concerned.  To  do  this,  it  is  sufficient  to  use  order  of  magnitude  estimates. 

* 

The  equation  of  motion  of  the  vapor  follows  from  eqs.  (2),  (8), 
and  (9)»  and  is  given  by* 

3y  , 

[git  2 ^ ^ * 2^  ‘ ^ 

The  velocity  in  the  vapor  is  certainly  smaller  in  magnitvide  than  the  velocity 
of  the  bubble  wall,  because  of  the  evaporation  which  takes  place  when  the 
bubble  grows  or  the  condensation  of  vapor  which  occurs  when  the  bubble 
collapses.  The  vapor  density  is  smaller  than  the  liquid  density  in  a ratio 
of  about  1:1000,  and  the  coefficient  of  viscosity  of  the  vapor  is  smaller  in 
a ratio  of  roughly  1:10,  The  pressure  gradient  in  the  vaptrmay  therefore 
be  assumed  smaller  than  the  gradient  in  the  liquid  by  at  least  an  order  of 
magiJ.tude.  An  estimate  of  the  pressure  gradient  in  the  liquid  may  be  made 

♦ ' ■ 

The  symbols  appearing  in  this  part  of  section  II  refer  to  the  vapor, 
unless  otherwise  indicated. 
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again  putting  the  liqtiid  velocity  at  the  bubble  wall  equal  to  the  bubble 
wall  velocity  ft  in  eq.  (17),  and  using  this  to  evaluate  eq.  (19)  at  the 
bubble  wall.  The  result  for  the  liquid  is  sln^jJiy  that 


!£iia 

6r 


r=R(t) 


o 

where  R is  the  radial  acceleration  of.  the  bubble  wall.  For  water  at  103  C, 

5/2 

the  maximam  radial  acceleration  of  the  bubble  wall  is  about  6x10  civsec 
(iit;.*  Ilg,  7),  giving  for  the  maximum  presstire  gradient  at  the  bubble  wall 
6 X 10^  d,ynes/cm  /cm,  or  about  .6  atm./cm.  The  maximum  pressure  gradient 
in  the  vapor  is  not  more  than  about  l/lO  of  this,  on  the  basis  of  the  above 

-3 

estimate,  and  it  occurs  when  the  bubble  radius  is  about  2 x 10  cm.  Thus, 
the  pressure  variation  in  the  bubble  is  at  most  of  the  order  of  10  ^ atm. 

But  the  pressure  itself  is  of  the  same  order  as  the  external  pressure  of 
1 atm.  It  is  thvis  clear  that  the  pressure  may  be  taken  as  uniform  within 
the  bubble, 

Pvap  = (26) 


For  order  of  magnitude  pm*po3es  it  is  sufficient  to  consider  the 
vapor  to  be  thermally  and  calorically  perfect: 


P = >»BT,  (27) 
e = c^T,  (28) 

where  B in  (27),  is  the  universal  gas  constant  divided  by  the  molecular 
weight  of  the  vapor,  and  in  (28)  is  the  specific  heat  (at  constant 
volume)  of  the  vapor.  The  heat  equation  for  the  vapor  becomes 


^°v^dt  + V • 7t1  = - p(7  . v) 

t • v)^  - 2v  • 7(7  ' v)^  (29) 

according  to  eqs,  (3)  and  (9),  if  the  radiant  heating  in  the  vapor  is  ignored. 
The  thermal  conductivity,  specific  heat  and  viscosity  coefficient  of  the  vapor 
are  about  an  order  of  magnitude  smaller  than  the  corresponding  quantities 
for  the  liquid.  On  making  the  same  approximations  for  the  heat  equation  (29) 
as  for  the  equation  of  motion,  we  obtain  an  approximate  relation 


- 18  - 


* 


which  reduces  to 

V^(kT  + p^5)  = 0 (30) 

if  the  velocity  potential  relation  for  the  velocity  and  the  uniformity  of  p 
and  k are  used.  Inasmuch  as  the  pressiire,  temperature  and  velocity  potential 
are  all  finite  at  the  origin,  the  only  solution  of  (30)  consistent  with 
spherical  symmetry  is  of  the  form 


kT  + = G(t), 


(31) 


The  velocity  potential  remains  undetermined  to  an  additive  function  of  time, 
which  can  be  so  chosen  that  C(t)  = 0 in  (31).  Eq.  (31)  then  yields  the 
relation 

0 = - ^ T,  (32) 

and  if  the  perfect  gas  law  (27)  is  used  in  (32)  it  gives  the  further  relaticn 


^ " JOB  * 


(33) 


Eq.  (33)  may  now  be  substituted  into  the  equation  of  continuity  (1),  giving 


= .y5v  = v .y>v0  = |v-^y, 
or 

J)B  — f • 

Eq.,  (34)  is  a diffusion-type  equation,  which  may  be  compar*ed  with 
the  heat  equation  written  for  stationary  vapor 


or 


if  = 


aT 

at* 


(35) 


The  function  in  eq.  (34)  replacing  the  thermal  diffusivity  D of  the  vapor 
in  eq,  (35)  is 


101 


k 


9 


(36) 
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If  the  -fraLper  aolsculss  nay  be  supposed  to  havei  say,  three  translational  and 
three  rotational  degrees  of  freedom,  then  =5  3B.  Hence  if  the  vJipor  were 
stationary,  ve  should  have  from  (35)  and  (36) 

\0\  z 3D,  (37) 

Eq.  (37)  will  be  of  the  right  order  of  magnitude  even  if  the  density  varies. 
Now,  the  thermal  diffusivity  D of  saturated  water  vapor  at  103  0 is  about 
.3  cm  /sec.  Hie  characteristic  difJ^sion  length  for  eq.  (34)  is  thus  about 
/4|j2(|t  js  2 \/t.  The  significant  time  for  the  103®  bubble  is  about  10  ^ sec, 
renighly  the  time  between  the  end  of  the  relaxation  period  (when  the  rate  of 
bubble  growth  becomes  significant)  and  the  time  of  the  velocity  maximum  (see 
Fig.  7),  giving  for  the  diffusion  length  i/4ljzfjt  a .02  cm,  av  about  six  times 
the  bubble  radius  at  the  velocity  maximum.  There  is,  therefore,  an  insignif- 
icant variation  of  In Ji  with  position  in  the  bubble,  so  that  the  vapor  density 
may  be  considered  a function  of  time  alone, 


Aap 

(38) 

Eqs,  (32)  and  (33)  then  show  that  also 

(39) 

«Up  = «'<*>• 

(40) 

It  is  not  legitimate,  of  course,  to  argue  fbom  eq.  (40)  that  the 
vapor  is  at  rest,  since  it  is  the  small  terms  in  ^ which  have  already  been 
neglected  in  arriving  at  (40)  that  determine  the  velocity.  This  difficulty 
dan  be  traced  to  the  normalization  chosen  for  However,  an  estimate  for 
the  velocity  is  readily  obtained  from  eq.  (38).  Consider  a sphere  of  radius 
r within  the  impor.  The  mass  of  vapor  in  such  a sphere,  by  (38),  is  simply 

ra(r,t)  = I irr^y?  (t).  (41) 

If,  now,  the  independent  variables  in  (4I)  are  chosen  to  be  m and  t, 
rather  than  r and  t,  eq.  (41)  may  be  witten  in  an  alternative  notation  as 

m = I n />(t)  r^(m,t), 


(42) 
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ths  cbaagd  correspoiidiag  to  the  adoption  of  lAgrangian  rather  than  Bulerlan 
coordinates.  Differentiation  of  (42)  with  respect  to  t»  holding  m fixed » 


shovs  ti»t 


or  since 


at 


is  now  the  vapor  velocity » that 


V 

vap 


(I'ft) 


^ SiS  . - . -X-  sLf 

3y>  at  ' 3/>  dt  • 


Inasimioh  as  the  origin  for  r in  the  above  development  was  arbitrary*  eq. 

(43)  inqjlies  a uniform  dilation  or  contraction  of  the  vapor  within  the  bubble, 
tMoh  is  consistent  with  the  uniformity  • Eq,  (43)  indicates  a radial 

velocity  in  the  vapor  at  the  bubble  vail  of  about  ,4  cm/ses  for  the  103° 
bubble  vdien  the  velocity  of  the  bubble  \iall  reaches  its  maxi  mum  of  about 
32  cm/sec. 

In  the  discussion  of  the  growing  vapor  bubble  in  a superheated  liquid, 
effects  related  to  the  velocity  of  the  vapor  will  be  neglected,  thus  incurring 
an  error  of  a few  per  cent  in  the  results.  In  the  discussion  of  the 
collapsing  bubble,  the  large  temperature  variations  at  the  bxibble  wall  may 
be  expected  to  eatiae  significant  changes  in  the  vapor  density,  and  so  the 
effects  of  the  vapor  veleoity  will  be  included. 

It  may  be  noted  in  passing  that  for  a uniform  (irrotational)  dilation 
€0*  contraction  of  the  vapor,  such  as  is  Indicated  by  eqs,  (38)  or  (43), 
the  viscosity  terms  in  the  heat  equation  (29)  vanish  identically. 


Boundary  Conditions, 

If  the  terms  involving  viscosity  are  omitted  Arom  the  basic  equations 
(1),  (2),  (3)  and  the  equations  of  state  are  substituted  in,  the  resulting 
aquations  are  cf  first  differential  order  in  v and  p,  and  second  order  in 
T,  In  order  to  match  solutions  for  the  liquid  and  vapor  across  the  bubble 
wall,  relations  must  therefore  be  provided  connecting  the  values  of  v and  p, 
and  the  value  end  normal  derivative  of  T across  the  interface.  Because  of 
the  assufiqption  of  spherical  symmetry,  no  further  relations  are  needed  If  the 
viscosity  teivas  are  retained. 
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Goaai«0r  a diff@rdatial  cone  F extending  from  xn©  cenxer  or 
the  vapor  bubble  (the  origin}  to  a point  indefinitely  far  away  in  tlw  liquid* 
idiose  generators  are  straight  lines  through  the  origin*  and  whose  differential 
orofis  section  at  the  bubble  wall  is  ^ • In  this  cone*  mass  is  conserved 


■at 


+ 


dS  Si  0. 


(44) 


Tbs  equation  of  motion  of  the  cone  is^ 

J ^2  dTi  a»^V2dS  = ^ n»PdS- 


i-p/  \j  -p 

and  ti^  energy  balance  equation  is 


^ a as, 


cr'dS  r + ^ n • /?z(q  + ^ ^ 


(45) 


J 


i dH  + 6 n « (P  • £ + m)  dS« 

r ^ r“ 


m 


In  these  equations)  £l  denotes  the  outward  drawn  normal  to  the  cone  in 
^ dS|  and  the  imit  nowaal  to  extending  away  freon  the  origin  in  ^ 

The  surface  integrals  on  the  left  side  of  the  equations  represent  transport 
terasj  these  must  be  added  because  the  elements  of  integration  do  not  move 
with  the  fluid*  To  evaluate  the  Integrals  we  shall  need  the  differential 
equations)  whiob  nay  be  written 


dS. 


af 

at 


+ V 


(47) 

(48) 

(49) 


Ms 

“g^  4 V • [y?V2  - Fj  = 0, 

[/  (e  4 I v^)1  4 ^ . [ y?2(e  4 i v^)  - P • 2 - W7T]  = 4f 
eq«  (49)  being  a oombinatlon  of  eqs.  (2)  and  (3)« 

Since  the  stress  tensor  does  not  account  for  the  molecular  forces 
responsible  for  surface  tension*  the  force  due  to  surface  tensions  must 
be  written  explicitly  into  eq«  (45))  and  the  corresponding  surfaos  energy 
into  eq.  (46),  For  a classical  treatment  of  the  theory  of  suyfae©  tension* 
sea  Joos*  Thsoratical  plroBics  (ifeffesr  FublishiE^  GOs*  Inc#*  New  York*  1934)* 
Chapt.  n*  seotion  8» 
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Tl®  equ&tl^s  (47)  i (43).  (49)  are  1*11  of  the  foriB 

^ t V . b = 0, 

a,  b aod  e being  suitable  scalars*  veot<»‘s  or  tensors.  The  integral 
relations  require  a knowledge  of 


(50) 


dt 


The  change  in  the  integral  | ad'D  during  a time  dt  is 


r -W+feit  L"* 


d l^T 


dt 


/. 


dt  - ] &’dt 
r<R 


r>R+Rdt  L 


a + ~ dt 

Ow 


“ I: 


ax'  - 1 adt;> 

'r>R 


(51) 


tdiere  a*  Is  the  value  of  a written  for  the  vapcarc  The  volume  element 
dx:>  is  here  considered  fixed  in  space.  To  first  (»rder  in  dt*  (51)  is 


d / adr  = / . (a*  - a)dX  + 

^ sj  T 


R<r<R+5dt 


dt[  f 

y<Il 


8a « 

at 


dt  + 


L 


ga 

at 


dt]. 


This  may  be  transformed  with  the  aid  of  the  relation  (50 ) and  the  divergence 
theorem*  to 

-fr  / adX  = / R(a»  - a)dS  + / o'dt  - 9 n • b»dS 

Qh  Jp  j P i/  P A.<r  “ 


L 


+ I Qd't 

r>R 


<p  n • bdS 
J r>ft 


= / R(a'  - a)dS  + i cdt  - <p  n • bdS  - f n '(b*-  b)d£ 

'^1  v/  p J p ^ 

= 2L  • [R(a’  - a)  - n •(b*  - b)]^  - ^ n e bdS  + J'  cdt  • 


(52) 
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Applying  (52)  to  sqs.  (44)“(A9)»  one  sAstains  in  ft  straightfcarsard  feshion 

[R(y!J»  -y?)  - a • (^‘2*  ~ 

[R( -y?v)  - a • (/>’2‘2*  •^22  “ “*  ■**  ~ “ * 

£SE^»(e +|v'2)  + lv2)l  I (53) 

- a *y<^‘2'(0‘  + 2 + 2 ~ 

+ P • 2 - k’^‘  + k^Tl  V . 

Beoai^e  of  the  spl^rioftl  8yimae‘toy« 

n . p a n[-  p + i (g  - J)1  (54) 

on  X (the  sane  relation  holding  for  P0«  ^ 'ftdft  (54)  in  (53)  leads 
to  the  results 

(R  - v)  =y5)  ‘ (R  - v‘)  f (55 ) 

P + ^ K p»  + (v  - v')y>'(R  - v»)  + 3^  (^  “ p)  **  3 

k ^ - k»  g-  = (R  " V’  )[L  + |(S  - V*  )^  - ^(R  - v)^ 

« 4 11  (M  . 8Z)  + ^ 2-  . Z)]  (57) 

3 W ar'  ^ 3/?  ^8r  r^^ 

at  r = Rf  \diere  in  (57)  ve  have  put 

L = e'  - e j (58) 

according  to  the  first  law  of  thernodynanics*  L is  the  latent  lieat  of 
evaporatiocx  at  the  bubble  vail.  Evidently;  eq.  (55)  esqiresses  the  can- 
servatloo  of  mooenttm  at  the  bubble  wall.  The  last  terns  in  (57)$  (58) 
represent  kinematic  (mass  transfer)  and  viscous  corrections  to  the  vapor 
pleasure  and  heat  transfer  relations  holding  at  the  bubble  vail, 

•»  mm  m 4 m 4^  MCft  ^ *1  *w  am  aX*  J Aw  s*  i«  ■■]  g 

A44KVAV  A9fp  4 i IT^  I I J | UajlO  UVUMAX  MXWIA  WX  «/UX  O W VAA  WXUMX  WJT  6»UX  wOO 

the  bubble  surftiee 

T»  s T.  (59) 

A toqperature  disoontiimitgr  would  iuq>ly  an  infinite  heat  eoDduotien  through 
the  Burfhoo* 
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The  momentua  condition  (3S ) , %Moh  laay  be  written 

V - Rii  - VI  - 7-ij» 

^ R 


(60) 


shows  that  the  liquid  velocity  at  the  bubble  wall  oan  dllTer  from  the 
bubble  wall  velocity  by  at  most  terns  of  relative  order  yo  « l/lOOO. 

For  all  practical  purposes  t (55)  nay  therefore  be  replaced  by 

V,.  (R)  - ft*  (61) 

llq 

As  has  been  neubloned  previously « the  vapor  velocity  Is  small  in  eomp^israt 
with  the  bubble  wall  velocity.  The  n<»ientum  transport  correction  In  eqj  (56) 
is  thus  approximately  ^’ft^*  For  this  to  coa^jare  with  *»  \diich  is  about 
1 atm  « 10^  e*g*s  units,  ft  must  equal  about  300  nv^seo,  i,d*  be  comparable 
with  sonic  velocity  in  the  vapor*  Such  large  velocities  will  not  be  con- 
sidered here,  so  that  the  momontua  tr‘anspo::'t  correotion  In  eq*  (56)  may  be 
neglected*  Similarly,  the  kinetic  energy  transport  corrections  to  the  heat 
transfer  relation  (57)  are  completely  negLlglble  in  conparlson  with  the 
latent  heat  of  evaporation,  (For  water,  L a 2 x 10^^  erg/gm.) 

The  viscosity  corrections  in  (56),  (57)  nay  be  evaluated  by  eqs* 

(17),  (43  )>  (60  )•  The  oonti'ibution  of  tlie  vapor  vanishes  identically* 

The  contribution  of  the  liquid  amounts  to 

in  eq»  (56),  and  this  divided  by  y?  in  eq*  (57),  The  net  effect  of  viscosity 
is  thus  to  increase  the  surface  tension  <7^  in  (56)  by  2:;^ft,  and  to  decrease 
the  latent  heat  L in  (57)  by  4^ft^R*  For  a vapor  bubble  grc»>ring  in 
water,  the  latter  oorreotion  is  entirely  negligible*  The  surface  tension 
increase  amounts  to  *2  and  *9  dyues/em  at  the  time  of  the  velocity  iiaximuffl 
of  the  103^  and  106"  bubbles,  respectively  (see  Figs*  7,9)*  This  is  of  the 
order  of  the  thermal  variation  in  idiloh  has  already  been  neglected.  For 
the  collapsing  bubble  § these  viscous  effects  become  i«^'Ortant  only  near  the 
point  of  collapse. 

With  the  ne^eet  of  the  klcematio  and  viscous  oorreotions,  the 
pressure  relation  (56)  reduces  to 


►'vap 


« + 2£. 
Pllq  ^ R • 


(62) 


Since  the  ten^rature  within  the  vapor  is  considered  to  be  unifccra,  the 
heat  transfer  relation  (57)  zeduoes  to 
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^ eq-  (43)*  this  inay  also  be  vrittsa 


When  the  vapor  is  at  rest  relative  to  the  interface » the  quantities 
P-vap  /^vap  equilibrium  vapor  pressure  and  density 

Pan  liquid  at  the  tesmerature  of  the  bubble  surface* 

But  ^»ien  a relative  velcolty  exists  they  viU  differ » Igr  an  amount  depending 
on  the  nature  of  the  liquid*  A relation  connecting  these  which  hftids  when 
the  vapor  may  be  considered  a perfect  gas  in  a state  of  cooqjlete  equilibrium 
^\an  “ ^®*^  given  by  Mathews on  the  basis  of  previous  work 

do) 

by  Plesset*  In  the  notation  used  herey  Mathews*  relation  is 


a -®  -- 

p /^eo  - . A * 
*^eq  c + R 


(64) 


vrtiere  o is  a characteristic  velocity*  related  to  the  velocity  of  sound  o 
in  the  vapor  and  the  specific  heat  ratio  Y of  the  vapor  by 


(65) 


The  naramater  a annAar^no  -in  Ha  mHoS  th«  MW  /«4 

• — *-*•-  — — — a r A A.  *l»  ViS^VAA  U 

of  the  liquid*  and  measures  the  fraction  of  the  surface  available  for 
evaporation  or  condensation*  For  non—polar  liquids  a is  near  unity  *but 
for  p(dar  liquids  with  bydr^en  bonding  a nay  be  appreciably  smaller* 

Pot  a water  surface  near  10®C,  o is  rsported  by  Wyllie^^^  ^ to  have  the 
value  *04*  This  es^erlnental  value  was  obtained  by  measuring  the  Mwg 
required  for  a san5)le  of  liquid  to  evaporate  into  a partial  vaciana*  Due 
to  ejqperinental  difficulties*  the  value  of  a for  water  has  not  been 
measured  at  higher  tenyieratures* 

If  the  value  o = *04  and  the  values  y = 1*33,  c - 5.  x 10^  om/eec 
are  used  in  (65 ) * they  give  for  the  characteristic  velocity  of  water 
0 « 7 i^Beo*  Correspondingly*  one  might  expect  a significant  deficiency 
in  vapor  pressure  and  density  to  occur  when  a vapor  bubble  is  growing 
in  water  at  only  a moderate  rate*  The  situation  in  the  case  of  a ooUapslng 
bubble  Is  even  more  orltical*  The  large  pressures  developed  within  the 
ooll^slag  bubble  beoanse  of  the  inaooesibility  of  bubble  surface  for 
oondensation  would  severely  limit  the  rate  of  odlapae*  and  mig^t  be  eaqeoted 
to  result  in  the  appearaaee  of  oondmsatloa  throughout  the  vapor* 
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These  oonoluslcmSf  however « do  not  appear  to  be  borne  ouit  by 
e^per^isesta  TbsiSf  the  pressure^Uinited  cciXspss  dosortbed  doss  not  ccour 
(radial  velooltlea  of  ooUapse  in  water  vhioh  are  certainly  in  excess  of 
25  n/seo  have  been  reported  by  and  condensation  heis  not  been 

observed  to  occur  within  collapsing  cavitation  bubbles « except  possibly 
near  the  point  of  collapse.  These  facts  indicate  that  the  value  of  c 
for  water  at  even  lO^C  may  be  aaieh  greater  than  and  possibly  point 
to  a significant  increase  of  a with  an  increase  of  teiiperatrare.'’^ 

He  shall,  therefore  assuaoia  the  velocity  o in  eq.  (64)  to  be 
BUffloiently  large  that  vapor  pressure  and  density  disorepanoies  may  be 
ignoredy  so  that  we  raj  take 

^vap  “ ^eq*  vap  *“  eq 

at  the  bubble  wall. 


n*oo  eqs.  (17)  y (61),  eq.  (21)  becomes 


O **  #2  ^ L.  *2 

R t 2RR*^  I BT 


(67) 


5*Mce 

the  bubble  wall 


by  sqs.  (62),  (66),  eq.  (67)  a*. 


HR 


S , I = Esai!L:£ 

* /'Hq 


2cr 


Aiq  ® 


(6S) 


an  equation  of  motion  for  the  radius  R of  the  bubble  wall.  Rq.  (68) 

(13) 

was  given  by  Plesset.'  refers  to  the  equilibrium 

vapor  pressure  of  the  liquid  at  ti^te  temperature  of  the  bubble  wall«1^^ 

5 ' ■■■I-  :-'  ' I :■  -T-  .--n  -t--:.  - l 

It  is  not  the  experimental  value  of  a which  is  questioned  hare,  but 
the  determiimtian  of  the  temperature  at  the  surface  of  evaporation.  In 
the  case  of  water,  the  rate  of  evaporatioi  is  large,  and  a steep  ten^rature 
gradient  develops  at  the  surface  (possibly  reaching  10^  or  10?  <’C/om). 


As  indicated  previously  in  tbs  discussion,  eq.  (68)  nay  be  otmsidered 
valid  so  long  as  h remains  small  in  oonparlscn  with  sonic  velocity- 
in  the  liquid.  A oorrection  to  the  equation  of  motion  (68)  %^oh  takes 
the  eosq^ressibiXity  of  thv  liquid  into  account  (up  to  terms  quadratic  in 
Voi|q)  has  been  ^ven  by  FyR.  Cdlmore  (I^ICIT  Report  26n4»  1952) 

on  iH  basis  of  tbs  Kitlcwo(^**6et]ie  feypothssis. 
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Couplad  id.th  tq*  (68)  is  ths  heat  equation  for  the  liquid » eq*  (24)* 

* 7Tl  = + 4,  (69) 

vith  the  boundary  oondition 

(70) 

(eqs*  ^63)j  (66))»  is  the  equilibrium  vapor  densii^ 

of  the  liquid  at  the  tamperatiEre  of  the  bubble  wall*  It  will  be  aaaudied 
that  Initially  the  teiaperature  in  the  liquid  is  uniform 

T(r,0)  * (71 ) 

together  with  the  initial  conditions  for  eq.  (68),  eqs.  (68)-(71 ) 
determine  the  problem  of  vapor  bubble  behavior. 


§ 


r=R(t) 
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III*  THE  IffiAT  PROBLEM 

The  system  of  equations  II  (69)  t (70 )»  (71)  whi.ch  define  the  heat 
transfer  problem  in  the  liquid  may  be  written 


5 [U  + V • VTl  - J 4(t)» 


81 

8r 


r=R(t) 


- F(t), 


> 


(1) 


jj 

where  D = — and 


T(r,0)  = T^, 


2 • 

R^  R 

V » 


eq 


J 


(2) 


At  r = w the  temperature  in  the  liquid  becomes  tiniform,  so  that  the  first 
and  last  of  eqs.  (1)  give  on  integration 


T^,=  T(oo,  t)  =T^+jq(t), 

if  q is  chosen  so  that  q(0)  = 0,  To  standardize  the  solution  it  is 
desirable  to  use  instead  of  T the  dependent  variable 


(3) 


e = T - T 


which  vanishes  at  r = « and  satisfies  the  system 


(4) 


+ vel. 


ae 

0r 


r=R 


= R'^F(t), 


> 


e(r,o)  = e{«> , t)  = 0. 


(5) 


- 29  - 


Convection  Solution. 

Because  of  the  boundary  condition  at  the  moving  bubble  vail,  it  is 
convenient  to  transform  (5)  from  Eulerian  to  Lagrangian  coordinates. 
Coordinates  appropriate  to  the  present  problem  are 


m s 1 [r^  - R3(t)3 , 


t = t. 


(6) 


The  Lagrangian  coordinate  m measures  essentially  the  mass  of  liquid  con- 
tained within  a sphere  of  radius  r about  the  center  of  the  bubble,  the 
liquid  density  having  been  assumed  constant.  In  terms  of  m and  t,  the 
system  (5)  becomes 


-S-  = 

0m  ^ D 0t  ^ 


4 _ 1 ^ 

” D at 

= F(t), 


ae 

0m 


|m=0 

©(m,0)  = 9(«  , t)  - 0, 


(7) 


These  eqtiations  can  be  put  in  more  tractable  form  by  introducing 
a temperature  potential  U,  defined  by 


e - M 
^ “ am  • 


(8) 


The  differential  equation 


I-  [r'i  e^LE  _ 1 22,  = 0 

am  ^ D at'  ^ 

may  be  then  integrated  once  with  respect  to  m 'to  yield 


AO  1 92  _ Tr+^ 

2 Dat"*'^^'* 


0m‘ 


where  J(t)  is  an  arbitrary  function  of  time.  From  eq.  (8), 

J'm 

© dm  + K(t), 
o 
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and  the  f\iriction  K{t)  may  be  chosen  so  that  J(t)  = 0,  and  also  so  that 
U(mtO)  = 0.  The  system  of  equations  to  be  solved  then  reduces  to 


L 6^  1 82  ^ 

8m2'“D8t- 


0, 


m=sO 


= F(t^, 


U(m,0) 


dm 


IZl^oo 


“ 0. 


> 


(9) 


The  diffusion  problem  thus  defined  can  be  solved  by  an  iterative 
procedure  if  the  assumption  is  made  that  U varies  appreciably  only  near 
m = 0.  This  is  equivalent  to  the  assTunption  that  the  temperature  variation 
in  the  liquid  is  localized  in  a thin  **thermal  boundary  layer'*  surrounding 
the  bubble  wall,  which  is  reasonable  if  (as  is  the  case  with  water)  the 
thermal  diffusivity  D of  the  liquid  is  sufficiently  small.  To  utilize 
the  assumption,  we  shall  rewrite  the  differential  equation  as 


4 2!k  _ 1 22  = (_4  . £d 
3^2  D at  R ) 5^2 


(10) 


and  treat  the  right  side  of  (10)  as  a perturbing  heat  source.  It  is 
convenient  to  use  in  (10)  a new  time  variable  ^ , defined  by 


'^  = J R^(t)  dt, 

in  terras  of  which  the  differential  equation  becomes 


a^u  1 au  ( ch 
8m2'5  W 


The  system  for  the  unperturbed  solution 


1 !!2  „ 


da" 


IffisO 

8U 


F(r), 

ra=co  *“ 


(11) 


(12) 


i 
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is  readily  solved  by  taking  the  Laplace  transform  of  the  system  with 
respect  to  the  variable  'i;;  , If 

u(m,s)  = / e“®^"  U (m,!:)  d "C  = ^ [U  J 

Jo®  ® 


and 


the  transformed  system  becomes 


f(s)  = 


d^u  s _ 


dm^ 


= f{s), 


dm 


ms=0 


= 0, 


with  solution 


Thus 


ia=oo 


rr 


j 


giving 


am* 


==^‘''  [^1  = f 'K  nlui  e’ 


f(s)  e 

-Vd 

d\ 

2 

= f(s) 

dm 

F(f  ) 

v/^ 

Elli. 

m 


4D(tf-  f) 


(13) 


(14) 


The  system  for  the  perturbation  correction  may  be  written 


a u.  . au. 

-ri  D8^  = 

om 


a^, 


dm' 


2 


(15) 


au, 


= D,  (a,0)  = ^ 


m=0 


= 0, 


m=«» 
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wnere 


G(m,r*)  = [ ^ - 1 


3^j 


2 2 
am'^  dm 


(16) 


If 


v(mis)  = g(m>s)  = ^[C*l» 


the  transfoimied  equations  become 


with  solution 


^-tv  = g(..3). 

on 


ijL 

'I 

dm*- 


m=0 


= 0. 


m=oo 


v=-|^  1 e 


'3  /§. 

-m  X 

6 


/ 


n ® 

® g(x,s)  dx  + e 


» -Xw  jv 

e ■ g(x,s)  dx 


■"•x/l 


r 


» "^v/l  1 

e ^ g(x,s)  dx  y t 

J 


so  that 


dm 


S . 


/s 


■) 

g dx  J. 


At  the  bubble  wall  m = 0,  this  reduces  to 


S2Z 

dm 


= - 1 e ^ ^ g(x,3)  dx, 

|m=0  «o 

giving  for  the  pertxirbation  temperature  at  the  bubble  wall* 


e,(o,r)  = - 


CO 


dx 


s/w  i 


It  is  assumed  in  the  following  dlsciiseion  that  all  limiting  procedures, 
changes  of  order  of  integration,  etc.,  are  permissible# 
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or  f^oia  (16) » 


a.(o,r)  =-“7^- 

V'iUD  > 


I 

./o 


dS  r 

- n^/2  J 


(^-J) 


o»  — 

xe 


8^U. 


a^, 


tac 


dx. 


(17) 


The  integral  equation  for  the  perturbation  potential  correction  vill 
be  given  for  cosq^letenesss  Its  I^iplaee  transfom  v(m;s)  can  be  written 
In  the  form 

/I 


T = -i./f  I 

‘ ' Jo 

£roxQ  which  one  obtains 

„ =.i  /I  r-AL. 
1 2J"  7o  /r^rj 


g(x,s)  dx, 


r 


, 4D(r  - J7  . ^ 4d('c  - 5 ) 


X 


/ / \ 

IS-' 


1 


lax= 


dx. 


The  solutions  given  above  in  eqs.  (t3)*'(17)  are  thus  far  exact, 
but  not  \28efUl  for  computation.  The  appraxlmations  to  be  made  depend 
upon  the  assumption  that  the  inflxience  of  the  heat  exchange  occurzd.ng  at 
the  bubble  wall  does  not  extend  far  into  the  liquid.  The  transform  of 
the  unperturbed  ten^erature  solution  also  be  written 


du  _ du 
dm  dm 


^ -'3. 
^JD 


in=0 


whieh  yields 


S^(m,7^ 


't'A  = -B- 


■/4sD 


pe^ojM^  -4d!^ 


(18) 


For  either  the  esq^anding  or  collapsing  bubble,  i6i^(0, 3 )!  increasing 

fUnotlon  of  3^  • Hence,  if  Q^iOtS)  W considered  negligible  for 
^ {<  "C),  one  obtains  from  eq.  (18)  the  inequall-ty 
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1 /_  -vx  f I ft  in  'y'  \ I f 

\iuy  u/j  ;;;  jv  tOf  '■'ti  w*iwi 


.a 


V \/4D('t  “ j 


^ V ^ 
\ u - 


n-  \ 

V 


Since  erfc(x)  <1/2  for  x > ,5#  the  oharaeteristio  diffusion  length 
la  terms  of  the  lagrangian  coordinates  may  be  taken  to  be  m = \/D(^  - 


for  the  unperturbed  teii^)erature  solution, 
perturbation  is 


The  parameter  of  the 


vhlch  in  the  region  of  significant  t^!5>erature  variation  is  therefore 
less  than 


Thus  the  perturbation  correction  may  be  expected  to  be  small  In  conqparison 
with  the  unperturbed  solution  vdien  3 \/D(  *^  - Xq)/R^  is  small  in 
ooD^rison  with  unitys  which  will  be  the  case  if  D is  sufficiently  small* 
The  perturbation  correction  must  vanish  when  0^  vanishes  (a.g*  for 
X,  < because  of  the  boundary  and  initial  conditions  on  the 

perturbation  equation. 

For  the  growing  bubble*  'ey  ~ < R^(t}(t  - t^)»  so  that 

3\/d(x  - X.J  3 /D(t“- tj 

'3  “ • 

R 

Tbs  significant  heat  transfer  for  the  103^  bubble  begins  at  about 

t^  = *15  millisoe.  The  bubble  raditis  is  2 x 10**^  cm  about  .06  ndllisee 
0 -3  2 
later.  Ttiklng  D“1«9x10'^  myaeOf  this  gives 

3 v/'SCt  - t J 
« .50 

near  the  beginning  of  the  103°  bubble  growth.  The  ratio  drops  asymptotically 
to  *20  at  later  times.  Thus  within  the  boundary  layer*  the  perturbation 
parameter  (rVn^)  “1  is  certainly  smaller  than  .50  d\u>ing  the  time  of 
signifieaat  heat  aarohange  at  tlu  bubble  wall*  For  larger  initial  8i;4)erheat8 


» 
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the  bubble  grows  fbstert  and  the  ratio  is  accordingly  soiM^^t  amallere 


For  the  collapsing  babble  considered  here*  the  ratio 


V 


is 


XDUch  smaller  than  unity  except  noar  the.  point  of  collapse* 

When  the  thin  ttiersial  boundary  layer  assumption  is  valid*  one  nay 
approxinate  the  perturbation  paraneter  by 


C' 


1 + 2a  ) . 1 „ m. 

R' } 


%ri.thin  the  region  of  significant  heat  transfer » and  to  a first  approocination 
neglect  the  term  (which  vanishes  in  any  ease  at  the  bubble  wall) 

in  eonq)ariaon  with  the  term  within  the  perturbation  solution 

integrals*  Thus*  the  perturbation  temperature  correction  at  the  bubble 
wall  eq*  (17)  beccoes  e^proxlinately 


9. 


I 


9 {0,^)=--^  / - 

% Vi 

or  by  eq.  (14)* 

" irD  R^(5)(X-f  )3/2 

/•i 

X f . 

-'o 


**  4D(f^  “»  ) Zx  . 

xe  ^ dx* 

o R-’  dx 


04> 


.3  4Dit  -f) 


X e 


dx 


Interchanging  the  order  of  the  last  two  integrations  gives 

io  (t-J 


eiCo.-S)  = -^  / 


R-^('ii-^)'''*-  Jo  Jo 


nD 


r 

*» 


= I 

" -4 


. a.w,rfa,*r.iy-go  • iSi-jL:; 

n^(r  -f  jj/2  ^ K-|  )‘ 


/ 

I > 


I 


(19) 
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nccOx'ding  bC  @qS  • (2  ) y (11 ) > 


F(r)  = i 


(20) 


is  proportional  to  the  rate  of  increase  of  the  nass  of  vapor  within  the 
bubble  I in  terms  of  the  time  variable 

^ j dt. 

F(i;)  is  therefore  negligible  tintil  a time  t = t^  or  “X  - when  the 

radius  of  the  bubble  begins  to  change  $ so  that  the  lower  limits  of 


X instead  of  0* 
o 


integration  in  eq«  (19)  may  he  taken  to  be  ^ ~ ^ 

For  the  collapsing  bubble*  |F(t.')j  Is  Itself  an  increasing 

function  of  'TU  « For  the  expanding  bubble*  R~Gi/t  as  t-*<»  (see 

Eq.  IV(62))*  so  that  after  an  initial  increase  P(o)  eventually  tends 
—I/2 

to  aero  as  '1^  ' j however*  for  a reasonable  choice  of  't  » the  product 


F(t) 


^ ^ in  this  case  becomes  an  increasing  function  of 
Hence*  the  perturbation  temperature  correction  given  by  eq.  (19)  is 
bounded  by 

t rr—TT  r f 


-e,(o,t)=f  / 

v'  V, 


n V>/ 


F(f)  . 

so  V 5 


— V 


R^(§) 


j/^ 


for  the  eaqpanding  bubble*  or  its  negative  for  the  collapsing  bubble 
last  Integtal  in  (21 ) can  be  transformed  bf  the  change  of  variables 

^o)(1  ~ x) 


(21) 

The 


1 


= + 


(0  < X < 1 ) 


3/ 
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giving  for  (21) 


r (lz\) 


3/2 


(20)  I this  ina7  b©  vritten  as 

’ /•^  >.v^  \3/< 


4DL 


3kR^(5) 


cr  if  u©  neglect  the  variation  of 


3k 


S<L 


lr(^)la'^ 

H^Ct)  • 

d J 


(22) 


d iii(^-  %) 


in  the  estimates  as 


4DL/J 


For  water,  the  factor 


y^eq 


(T)  in  eomparison  with  that  of 


(r^  ) 

o ' o' 


W 


d In  R^(^) 
d ln(3  - ^ ) 


(23) 


4DL/> 


£a 


3k 

in  (23)  has  the  value  ,066°C  at  T = 50®C,  and  r50®C  at  T = 100®C. 

In  the  eas©  of  the  ooUapsing  bubble,  the  ratio  of  logarithmic 
derivatives  in  (23)  is  small  except  near  the  point  of  collapse,  when  it 
increases  rapidly  in  absolute  value*  Since 


± 


(t-  X^)- 


for  crude  bound  for  j©^l  for  the  collapsing  bubble  is 


4DL^ 

|G,|  <— ^ 


3k 
4CI.  /> 


A d InC^-  -CJ 


'V 

O 


3k 


3BL 


( 

In  -! 


3 

Xs. 

^ Y? 

V J 


d In  yfi) 


d ln(^  ~ ^ 


(24) 


This  haa  about  the  wlue  *07^0  x 6.9  ts  »5^C  £nr  the  bubble  considered  here 
vhen  the  bubble  radius  has  dropped  to  1/10  of  Its  initial  'Talue*  which  is 
ftother  than  the  collapse  was  actually  followed  (see  Fig.  11  )•  Howe'ver* 
the  actual  temperature  rise  (see  fig.  12}  at  the  end  of  the  tine  period 
considered  was  40°C»  in  cooparlson  with  which  the  perturbation  oor3reoti(»i 
estimate  is  eon^letely  negligible. 

For  the  expanding  bubble*  the  ratio  of  logarithmic  derivatives 
appearing  in  sq,  (23)  vanishes  at  rises  to  a maximum  cf  less 

than  *75  at  some  value  of  i then  drops  asyn^jtotlcally  to  the  value 


1/2.  The  decrease  after  the  maxlsum  is  sufficiently  gradual  that  it  can 
be  d^idnated  by  a factor  ('^-  ® small  value  of  c»  Thus*  a 

bound  on  |e^  | for  the  expanding  bubble  la  given  hy 

/V  \l/2*'C 


, I f V 

V - 3k 


uc 

o 


372 


4^, 

3k 


3/2  - c 


d In 
d ln('!^- 

d In  R^(y) 
d ln('^ 


X 


r 

o 


d in(^  - 

(/-  d'5 

Vz 


(25) 


Taking  for  the  103®  bubble  to  correspond  to  t^  = .15  milliseo  (see 
F5.g.  7)»  we  find  for  the  ratio  of  logarithmic  derivatives  in  (25)  a 
maxiinuiii  of  *71  at  t ~ .19  milliseo**  ond  the  value  c » .IS.  The  bound 
on  |9^  { thus  has  a maximum  of  about  •26®G  at  .19  milliseo*  and  drops 
asys^totically  to  .18®C.  For  the  106®  bubble  (Fig.  9)  the  choice  t^  = 28 
sdcroseo*  gives  a maxlimim  ratio  of  *74  at  34  mlcroseo.  and  the  value 
c » .25.  Since  4DLy6'^^/3k  s *60®C  at  106®Cf  the  bound  on  |©^|  for  the 
106®  bubble  has  a maxlinuia  value  •36®C  and  an  asyn^totie  value  .24®C. 

The  conclusions  dravn  £rm  the  analysis  above  of  the  pertwhation 
parameter*  and  f^om  the  bounds  derived  for  the  pertiurbation  teo^erature 
correction  at  the  bubble  wll*  are  thus  in  agreement.  They  indicate  that 
the  thiji  thermal  boundary  layer  aeswpti^  is  valid  when  the  vapor  bubble 
grows  or  collapses  rapidly*  and  also  my  be  considered  valid  when  the 
bubble  is  nearly  in  eqtailibriimi.  The  assun^tioo  beoooes  somewhat  critical* 
in  the  case  of  water  above  its  boiling  point*  vdien  the  bubble  growth  occurs 


rats  ^ 2uffieisnt2y  rapid  that  significant  haat  transfer 
occura  at  the  bubble  wall^  but  slow  enough  for  the  liquid  to  partially 
adjust  to  this  transfer.  The  error  incurred  at  this  critical  stage* 
however*  does  not  amount  to  more  than  a few  tenths  of  a degree*  and  remains 
a small  fraction  of  the  actual  temperature*  It  is  sufficient*  therefore* 
to  use  the  unperturbed  temperature  solution  in  the  dynamic  problem*  at 
least  for  the  case  of  water* 

In  terms  of  the  original  time  and  radius  variables,  the  unperturbed 
temperature  solution  (13)  is 


®Q(r»t)  T(r. 


Jo  /FTI 


r=^(x) 


dx 


R^(y)  dy 


X 


exp 


^ 36d  I R^(7)  dy 

y •St  ^ 


which  reduces  at  the  bubble  wall  to 


©o(K»t)  = - 


The  ten?)erature  gradient  at  the  bubble  wall  is  given  in  terms  of  the 
e^pnration  rate  by  eqs*  (1J*  (2)*  so  that  for  the  present  problem*  the 
unperturbed  solution  beoomes 


(26) 


PlffUBioa  Solutions* 

It  is  of  interest  to  ooo^re  with  the  convection  solution 
pres^tsd  above*  a diffusion  a^roodsiatlOT  (for  the  identical  physical 
problem)  proposed  hy  Forster  and  Zuber.^*^*^^^  This  apprcKimation  may 
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bs  dsTslepsd  in  the  foUswing  Bsarjisrs* 

The  heat  equation  (5}  for  an  Inecanpreasible  fluid  nay  be  vo'ltten 


H = ^ • [k^Q  - V /c^Ol . (27) 

In  this  form#  one  may  identic  the  last  term  on  the  right  as  representing 
the  rate  of  heat  influx  into  an  elementary  volume  fixed  in  space  due  t& 
transport  the  fluid#  If  the  fluid  velocity  is  sufficiently  small# 
the  transport  effect  may  be  neglected  in  eamparlaon  idth  conduction  effects# 
represented  by  the  first  term  on  the  right.  With  this  neglect#  the  equation 
bee^oes  a diffusion  equation# 


^ ® D at* 


®v 


(28) 


The  initial  and  exterral  boundary  conditions  f^  eq#  (28) 

©(r,0)  = e(o«  # t)  0 (29) 

are  the  same  as  for  the  convection  equation#  The  boundary  condition  at 
the  moving  vapor  bubble  wail 


- k • 4nR 


0r 


= in  (| 

^ dt  3 / eq 


(30) 


A — *•  \ V / 


is  also  the  same  # but  becomes  extremely  difficult  to  apply  in  £!ulerlan 
coordinates#  and  seme  sort  of  physical  or  mathematical  artifiee  must  be 
resorted  to  if  a solution  in  closed  form  is  to  be  obtained#  The  approach 
•^o  the  problem  given  by  Forster  and  Zuber  consists  of  treating  the  heat 
''Exchange  at  the  bubble  vail  as  though  it  were  due  to  a moving  spherical 
heat  source  (for  the  collapsing  bubble#  or  a heat  sink  for  the  expanding 
bubble)  in  a stationary  infinite  fluid# 


The  presentation  given  here  is  not  that  of  the  authws#  Tlie 
original  presentation  of  Forster  and  Zuber  (ref#  (7))  la  quite  brief | 
the  second  (wf#  (14))  by  Forster  purpsrts  to  give  a more  detailed 
treatment  of  the  problem#  but  aottially  treats  a different  problem# 
Neither  paper  gives  an  adequate  analysis# 
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They  start  with  the  elementary  solutlc^* 


© = 


4JTrr>  v/^5t 


jfaii: 

4Dt 


4Dt 


(31) 


of  equaticais  {2B)t  (29)  for  the  temperatu-re  difference  0 in  the  fluid 
at  any  radius  r and  any  time  t > 0 due  to  an  instantaneous  spherical 
heat  source  at  r = r’»  t = O.  The  total  heat  liberated  is  2^c^Q* 

This  is  to  be  related  to  the  quantity  of  heat  transferred  out  of  the 
bubble  by  condensation  during  a time  interval  dx  vdiile  the  bubble  radius 
is  r’  = R(x).  The  actual  heat  h transferred  cut  of  the  bubble  at 
time  t X is  given  by  the  right  side  of  eq*  (30). 


h = “L  ^ v R^(x)/»gq1  dx,  (32) 


and  this  is  also  the  heat  transferred  into  the  liquid  at  t = x.  ]ji 
accordance  with  the  Forster  and  Ztiber  assumption  of  an  Infinite  medium^ 
however;  the  elementary  heat  source  associated  with  the  solution  (31 ) 
releases  its  heat  not  only  to  the  fluid  outside  the  shell  r = rS  but 
also  to  that  part  of  the  fliiid  inside  the  shell.  At  th©  instant  of  release ; 
half  of  the  heat  appears  inside  the  shell.  Therefore;  if  the  solution 
(31 ) is  to  correspond  to  the  heat  release  outside  the  shell  at  the  moment 
of  liberation;  the  h of  eq.  (32)  must  be  equated  to  only  half  of  the 
total  heat  liberated;  giving 


This  choice  iutroduces  on  error  at  later  times  in  the  final 
solutlcm  for  two  reasons:  First  of  £lL1;  the  beat  flowing  through  a later 
shell  is  no  longer  Just  tteit  due  to  condensation  at  that  time;  but  still 
has  a oontribution  frcm  the  heat  diffusing  outward  from  the  shell  r = r’. 
This  relaxation  effect  la  minimal  if  the  radial  velocity  of  the  bubble 
wall  is  large;  and/or  the  thermal  conductivity  of  the  (stationazy)  fluid 

Eq*  (31)  is  readily  obtained  by  the  operational  methods  already 
used.  It  dlixers  from  the  standard  solution  given  by  Ci^slaw  and 
Jaeger.  Conduotl<m  of  Heat  in  Solids  (Oxford  Univ.  Press#  1947);  p#  219? 
by  a factor  of  tv©  because  of  ilia  choice  2 ^ c Q;  rather  than  yo  c Q* 
for  the  total  heat  liberated.  Tk»  transform  of^(31 ) is  given  below.^ 
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/•w..,.. 

... 


is  isa&lle  Ihif ; these  sre  not  the  conditions  of  of  the 

diffosian  modsl*  Saoondllya  ths;^  is  et.ill  th^  contvibvition  tu  the  heat 
extent  of  the  liquid  at  later  timaa  by  the  hsat  liberated  vithln  the 
shell  r - r^«  heat  released  within  the  shell  ultiaately  reverses 
iti  ^reetion  of  flow  and  adds  to  the  heat  ocaitent  witside  the  shells 
This  effeot  can  be  eliiainated  by  a aisqtle  davioe  to  be  described  belcfw« 
EUisimticn  of  the  relaxation  effect  is  a more  Involved  matter. 

Bie  Forster  and  ZOber  solution  is  obtained  by  replasing  t fey 
t - X in  eq«  (31)  (eorresponding  to  a heat  release  at  t = x»  rather 
than  t = 0)  and  r'  by  R(x)»  using  (33)  to  define  Q»  and  integratiag 
the  result  ever  all  sources  from  x ” 0 to  x=tg  This  gives  for  the 
tsj^srature  dlfferenoe  at  the  bubble  wall  r “ R(t)f 


®(R,t)  = - 


— . 6,  f*  ^ ^ 

3/0  Wo  R(t)  R(x}  \^-x 


mu  - 

4D(t-a) 


[■»(*).♦  H' 

4E(t-*j 


• t 

J ‘ 


Memlag  thaty  for  the  dynawi®  problfin  under  ocnasideratlon  (the  growing 
bid3ble)y  the  second  e:q>onent^al  fa  (34)  can  be  neglected  and  tjie  first 
by  t^tyg  Forster  and  i|ttber  further  reduce  the  above,  equation 

tc  "-'..  f' 

/'"  ./  ■ ,y  ■ " ' ■ ' ' ‘i  :: 

A.  t Ju  I fCtf  ' Afl  - - _ 4. 


; . ^Ryt)  - lu_  . (35) 

3y>o^  \/«B  R(t)  R(x)  V^ 

r A soaew^^t  different  ap^roash  to  the  diffusion  problem  with 
moving  V/bimdaryy  for  idiioh  the  difiloulty  conneoted  id.tb  the  heat 
release  iheide  the  boundary  does  not  arise y is  the  fallowings  Lat  it 
be  e^P'^ed  that  at  t.  r Oy  a (mathematiool)  shutter  cjiens  at  r " r®y 
sug«  as  to  |m#g  thri^h  the  shell  a presaribed  quajitity  of  heat  hy  than 
instantaneously  ibloeee  Itgain.  Ihe  effect  ©f  the  moving  bouxidary  can  then 
olr^iMdf  *§  above  y by  Inte^ating  avsr  successive  shutters  at 

"r’  , S?  R|^.)y-  '0  < X < "te 

. .v  -Bo^  ..pa^irs  Fors'^r  attempt  to  justify  tiw  lorsMle.  (35)  as  an 
to  the  soiUtioh  of  the  diffusion  pr^lemy  and  as  an 
sigsrc^i^te  sointldu  t®  the  osn-e^ctlcsn  problwsa.  Bbe  ocavsotl.^,  sppp^- 
aatien  availabley  andy  in  fact?  referred  to  is  both 

. lianers. 
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The  ’’Bhu'fe't'SJ*  ss.21  bs  obbainscl  by  Tiot-ing  ■ths  hes't 

flow  through  a surface  element  boundirkg  the  fluid  is  determined  (as  for 
the  oonYeotion  solution}  by  the  ten^rature  gradient  there*  The 
appropriate  boundary  condition  for  a shutter  at  r-r',t*0  will  thus 
be  of  the  form 


ae 

dr 


= Q‘<?  (t). 


(36) 


r=5*’ 


The  constant  Q*  can  be  related  to  the  beat  h passed  by  the  shutter 
by  a determination  of  the  total  heat  outside  the  shell  at  a later  time 
(t  > 0). 

The  system  of  equations  to  he  solved  is 


= 1 SS 
® D at 


0(r,O)  = 0(«>  y t)  = Of 


96 

8r 


rsqr* 


(r  > r’ )» 


= "Q‘<^  (t). 


(37) 


% putting  w * ^ [©1  and  taking  the  laplace  transform  of  the  system 
with  respect  to  tf  one  can  reduce  it  to  an  equivalent  system 


(rw)  = I (rw), 
dr 


w(«»  , s)  = 0» 


& 

dr 


r=r' 


with  solution 


• Qt 


w(r,s)  = — • i e 


. 1 


(r  > r»). 


The  heat  liberated  outside  the  shell  r « r*  is  given  by 

.2  _ ,_„,2  ^ ^ D 


A.  “ J ^ ^ ” 4nr’“‘/£?  o^Q*  • J » 


(38) 
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i»ei 


h = C„0*D, 


80  that  (38)  becomes 


w(r,s)  s: 


1 e 


4y>o^D  r 


1 + r« 


(39) 


vhich  has  the  inverse  tranafom 


Q * 




4Wy«e^rr’  v/^rDt  ^ 


e 


p-r 
■ 4Dt 


n2 


-fit  _ (£-  . 1 ) 

,2 


erfo 


r «♦  r’ 

V® 


\ 


(40) 


^he  analogue  of  the  moving  source  solutlcm  (34)  can  be  obtained 
fran  (40)  by  again  writing  t - x for  t»  R(x)  for  r*,  R(t)  for  r» 
using  eq»  (32)  for  h»  and  integrating  over  x«  The  moving  shutter 
analogue  is  thus 


© = - 


— f 

_ x/tD  '^C 


/3  ) ife 

dx  ^ y^ea^  ^ 


4D(t“x) 


3yoe^  x/irD  '^o  R(x)  R(t)  s/t-i 


i 


e 


%ai  + sM_i 

r5)  ® 


erfo 


^ filthSlad 


(41) 


This  solution  has*  of  course » the  same  difficulty  with  the  relaxation 
heat  flow  that  the  moving  sou.ree  solution  hasy  but  the  problem  associated 
with  the  heat  flow  from  within  the  moving  boundary  baa  been  eliminated* 

It  Is  posslblS)  in  principle » to  eliminate  the  zelaxatlon  heat 
£l.ov  from  the  shutter  solution  by  using  the  correct  boundary  condition 
at  the  bubble  vall«  This  can  be  done#  for  instance » by  leaving  h 
undetermined  in  (40)  • but  summing  the  gradients  of  the  elementary  shutter 
eolutlone  over  the  variable  x to  obtain  the  ten^ratur#  gradient  of 
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the  fina3.  golutiwi  la  teruta  of  the  unknovm  differential  h(x)«  ^Jr  using 
eq,  (30)  to  specify  the  temperature  gradient  at  the  bubble  wall,  one 
obtains  an  integral  equation  to  be  solved  for  h(x}s  The  same  procedure 
can  be  carried  out  for  the  moving  source  solution,  leading  again  to  an 
integral  equation  for  h(x)«  It  is  apparent,  however,  that  a more 
appropriate  procedure  would  be  to  return  to  the  original  equations  and 
to  attempt  to  solve  them,  using  the  ecrreet  boundary  conditions  from  the 
beginning* 


The  moving  source  solution  (34)  and  the  shutter  analtgue  (4I ) 
differ  fb(a&  one  another  in  the  last  terms  of  the  respective  integrands* 
Since  this  difference  accounts  for  the  false  heat  flow  fp<ai  within  the 
Interior  of  the  bubble  in  the  case  of  the  moving  source  eolution,  it  may 
bs  expected  to  become  iiqjortant  wJ^never  relaxation  effects  become 
Important,  i«8*,  when  the  thermal  diffusivity  of  the  fluid  is  large*  This 
may  be  verified  indlrectlly’  by  showing  that  the  solutions  become  Identical 
when  the  diffusivity  is  small*  (A  direct  verification  will  be  given 
below  for  the  case  of  the  growing  vapor  bubble*)  Perhaps  the  easiest  way 
to  show  this  is  to  examine  the  Laplace  transforms  of  the  respective 
elementary  solutions.  If  the  transform  of  (40)  is  denoted  by  w,  as 
above,  and  that  of  (31;  by  u,  one  has 


(r  > r’), 


u(r,s)  = 




K^/>  o„D 


1 

r 


-{rtr* 


(r  > 0)* 


In  the  equation  for  u,  the  substitution  Q = h/^  c has  been  made. 

It  will  be  Observed  that  for  r > r*,  r*  /»  » 1,  both  reduce  to  the  form 

V D 


vbich  h&8  th8  inverse  transform 


© = 


rr*  \/iiDt 


e 


■ ir-rVi* 

4Dt 


contribution  to  the  boundary  tesqierature  by  tld.8  elementary  solution  is 


e 


JL 


) dx 


R(x)  R(t)  v'^t^c 


(42) 


i*««  the  leading  term  of  the  solutions  (34) » (41) • 

The  solution  (42)  is  a valid  a^prcadmatiou  to  (34)  &nd  (41) » 
however  f only  for  small  D«  For  large  C»  the  reznaining  terms  in  (34); 

(41 ) become  important}  as  C tends  to  finity;  the  integrands  of  both 
(34)  sad  (41)  tend  to  zeroi  in  fsot*  for  x ^ t«  £q»  (42)f>  on  the  other 
hand;  tends  to  the  Forster  and  Zuber  approximation  (35)  with  large  Ds 
It  is  to  be  concluded  from  this  that  the  Forster  and  Suber  approximation 
is  never  a valid  approximation  to  either  of  the  diffusion  solutions 
(except  in  the  case  of  a quasi-stationary  bubble  R(t)  » constant).  IJhen 
it  is  permissible  to  neglect  the  second  ejqzonentlal  in  eq.  (34)  in 
comparison  with  the  first;  the  first  eaqjonential  cannot  be  set  equal  to 
unity}  and  when  the  first  e^qponential  epproaohes  unity;  so  does  the  second. 


General  Ccmparlaon  of  the  Convection  and  DLffusleati  Solutions. 

The  temperature  solutione  presented  abov($  have  not  been  restricted 
thus  far  by  a radius-time  relation;  and  so  may  be  compared  for  ar^ 
assumed  behavi.or  of  the  boundary  consistent  with  the  heat  problem  (such 
as  may  be  achieved;  for  instance;  bF  keeping  the  teiiperature  of  the  bulk 
liquid  at  but  varying  the  extsrml  pressure). 

For  a growing  bubble;  the  convection  solution  (26)  predicts  a 
teoperature  drop  at  the  bubble  wall  which  varies  inversely  as  the  square 
root  of  the  thermal  diffUsivity  D of  the  liquid,  other  factors  (suoh 
as  the  specific  R(t)  behavior)  being  held  constant.  The  diffusion 
solutions  (34);  (41)  predict  a smaller  drop  than  the  convection  solution 
for  all  D.  The  discrepancy  becomes  most  marked  for  small  D;  vhen  the 
oonvectiom  drop  beoomed  large  but  the  difibsion  drops  tend  (depending 
vpm  the  lav  of  growth)  to  vanish. 
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7hese  predictlmts  are  qualitatively  xmderstandable  on  the  baais 
of  the  physical  models  involved*  In  the  case  of  the  oonveotion  model* 
the  heat  source  (or  sink)  la  alvaya  located  at  the  seune  fluid  elements* 
those  at  the  bubble  vail*  l%us*  a decrease  in  the  diffusivity  has  the 
effect  of  eoneentratixig  the  region  vithin  which  algnifioant  heat  transfer 
occurs  neeurer  to  the  boundary*  and  correspondingly  the  ten^erature  drop 
will  be  greater  there  (if  the  bubble  is  growing*  or  the  temperature 
rise  greater  there  \dien  the  bubble  is  collapsing)  for  small  D than  for 
large  D*  In  the  case  of  the  diffusion  model*  on  the  other  hand*  the 
fluid  remains  stationary  while  the  beat  soiu'ce  sweeps  through  it*  a 
decrease  in  the  diffusivlty  here  may  ultimately  be  esqpeoted  to  have  the 
effect  of  insulating  the  successive  elementary  sources  tram  one  another* 
such  as  to  prevent  any  acoumulation  of  heat  from  taking  place* 

l/hen  the  dlffusivily  of  the  fluid  is  large*  the  relative  inqpert- 
auoe  of  Qonveotitm  effects  in  comparison  with  diffusion  effects  should 
diminish*  As  has  been  pointed  out  above*  the  discrepancy  between  the 
convection  and  diffusion  solutions  is  most  marked  at  low,  rather  than 
high  'mlues  of  the  diffUslvity*  However,  both  the  omveetion  and  diffusion 
solutions  presented  hero  cease  to  be  valid  when  the  diffUalvlty  is  too 
large*  so  that  a caparison  in  the  limit  of  large  D is  not  sj^aningful* 

Conmarieoa  for  gree  Bubble  Growth* 

When  the  bubble  growth  is  not  forced  by  exterital  pressure 
variations*  it  beccews  limited  eventually  by  the  heat  transfer  at  the 
bubble  wall*  The  physical  relations  bolding  In  the  heat-limited  growth 
will  be  discussed  later,*  but  may  be  briefly  related  here*  The 
evaporation  at  the  bubble  wall  necessary  for  bubble  growth  forcjes  the 
tea^rature  of  the  licjuid  thera  down  toward  the  boiling  point  of  the  liquid 
at  the  external  pressure*  If  the  boiling  point  is  denoted  by  5^*  and 
the  temperature  of  the  liquid  by  T^*  the  late  growth  of  the 
bubble  must  then  be  such  as  to  satisfy  the  asymptotic  relation 

e - - (T  - T.  )*  (43) 

o 0 


See  the  discussion  ct  the  aiQmgitotio  phase  of  bubble  grow^  in 
aeettoa.  IV* 
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For  every  tenperature  solution  presented  here*  the  relation  (43)  restricts 
the  late  bubble  growth  to  a lav  of  the  form 

R{t)  as  t • (44) 


The  rate  of  bubble  growth  is  then  essentially  specified  by  the  value  of  C« 
Let  the  parameters  S»  A be  denoted  bF 


(T„  - Ifc) 


Then  S is  related  to  A by  an  equation  of  the  form 


(45) 


S = I(  A),  (46) 

%diere  I(A  ) denotes  the  temperature  Integral  involved*  For  the  oon" 
veetlon  solution  (26) » 1(A)  is  given  by 

i(x)  =Av/5  P — = — r~A>  (47») 


for  the  moving  source  solution  (34)  by* 

fZ 


r 


e 


hv5  _ 


. 1-v^j 


= 4X* 


^1  -A  , erfo(A)^  » 

(1  - A \/ff  + •••)  as  A 0 (D  » 

2\/n  (1  - + *.*j  gg  A ”*  w (D  •*  0)| 

2A 


for  the  moving  shutter  solution  (41 ) by 


UTj>) 


« 


The  integrals  (47b )»  (47c)  are  evaluated  in  Appendix  B* 
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iW  =x  j' 


J 

\ 


( 


9 


-i«y? 

J 

as  X 0 “*)» 

as  ^ 00  (D  0)j  C4'Q) 


2)?  (1  - — ) 


2\/n  (1  - “%+  •••) 

A* 


and  for  the  Forster  and  Zuber  apprcsdoation  (35)  by 


ia)=>  = (47d) 

s/T^X 

It  vill  be  observed  Srasn.  (47b )»  (47e)  that  tfcie  moving  souroe 
soXuticn  and  shutter  solution  become  asyn^totlcally  equal  in  the  region 
of  large  X (ios*  of  small  D»  for  a given  valve  of  C)»  but  differ  by 
a factor  of  two  uhen  A is  small(D  large).  This  factor  of  two  is  to 
be  attributed  to  the  ccnti'ibutiou  to  the  heat  content  of  the  liquid  by  the 
false  heat  flov  from  within  the  bubble  surface » iMch  occurs  in  the  case 
of  the  moving  source  solution  (47b}*  This  just  doubles  the  e]q>eoted 
content  (and  temperature  difference)  for  large  D.*  It  may  ctlso  be  noted 
that  the  Forster  and  Zuber  approocimation  (47d)  behaves  differently  in 
all  ranges  of  A moving  souroe  solution  (47b)  which  it  is  suppesed 

to  8q>proocimats* 

Ihysicallyf  lubble  growth  with  a given  value  of  C but  various 
values  of  D can  be  obtained  by  choosing  llq;id.ds  with  differing  thermal 
diflUsivlty«  and  adjusting  the  si;q>erheat  in  each  to  give  the  specified 
rate  of  growth*  If  one  concentrates  on  a given  liquid  (with  fixed  D) » the 
parameter  which  varies  with  the  superheat  becomes  G;  so  that  the  analysis 
of  the  varlaus  tes^rature  soluticuis  for  a given  liquid  is  to  be  made  on 


The  actual  heat  ocntant  here  is  negative » corresponding  to  the 
heat  less  from  the  liqpald  at  the  surface  of  the  growing  bubble  due  to 
evaporation* 
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tile  basis  of  the  clsnaemio  situation  vhiah  occurs  in  tiie  raspeotive  idigrsleal 
models  I aad  its  dependence  on  the  superheat* 

When  the  lav  of  ^ovt^h  (44) » applies*  the  tesqperature  at  the 
bubble  vail  has  dropped  practically  to  the  boiling  texqierature  of  the  liquid 
at  the  external  pressure*  Tbs  bubble  growth  is  maintained  by  a differential 
tenrperature  (and  therefore  pressure)  effect  which  vanishes  with  the  radial 
velocity  of  the  babble  vail*  and  which  is  negligible  in  ccoiparison  with  the 
t^erature  difference  - T^^.*  The  physical  constants  o^»  ^*/^eq 
and  D oharaeteristio  of  the  liquid  nay  here  be  given  their  values  at  the 
boiling  point  The  asy^totlo  rate  of  bubble  growth  (the  constant  C 

in  eq*  (44))  for  a given  model  is  then  determined  by  the  relations  (45)*  (46) 
and  the  superheat  T^  - T^  of  the  bulk  liquid* 

. The  thermal  relaaoitlon  effects  which  make  the  boundary  condition 
at  the  Ixibble  wall  inaccurate  * in  the  case  of  the  diffusien  soluti^s*  become 
important  if  tiie  bubble  wall  moves  too  slowly*  SincB  l(X)  is  in  all. 
oases  an  inoreasing  funotlon  of  \ * this  means  (aoeording  to  eq«  (46)) 
that  the  diffusion  solutions  do  not  represent  tiie  diffusimi  model  at  low 
superheats*  At  larger  superheats#  those  solutions  beoome  adequate  representa- 
tions of  the  diffusion  model*  But  since  tiie  radial  velocity  of  the  bubble 
wall  increases  with  the  superheat#  the  diffusion  mcdsl  itself  becomes  non- 
pbysloal  at  larger  superheats* 

The  adequacy  of  the  eonveotion  solution  may  be  determined  ft>om 
eq*  (^)s  Viten  the  sayiwjtotio  law  of  bubble  growth  (44)  applies#  the  ratio 
of  log^ithmio  derivatives  appearing  in  (25)  is  equal  to  1/2*  !nie  perturbation 
oorrecticn  to  the  ocnvectica  solution  (26)  is  therefore  not  lar^r  ti^n 
about  *2*^0  for  water#  and  aooordlngly  is  negligible  for  all  but  the  lowest 
siqMTheats#  once  the  vapor  bubble  grovtii  has  reached  the  asynptotio  stage* 

One  may  therefore  oonslder  the  oonvect<ion  s elution  to  be  accurate  in  this 
phase  of  bubble  growth. 

A plot  of  the  values  of  C predicted  by  the  convection  solution 
(47a)  and  the  moving  statter  soluticn  (47e)  for  varying  degrees  of  superheat 
'•All  be  found  In  Fig*  2#  fw  the  ease  of  water  at  one  atmosphere  external 


If  the  differential  temperature  effect  were  ignored#  the  bubble 
growth  rsquircid  by  the  theory  would  appear  paradoxioal* 
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pressure*  ^le  paraiseters  /O  » e . L»  />  and  D have  all  been  given 
their  values  at  = ICO  C*  The  ten^erature  intep'al  (47e)  fer  the 
aoving  shatter  solution  vas  integrated  nunerloall7B  The  breakdown  of  the 
diffusion  model  for  water  is  clearly  hhown  in  Eig*  2*  The  diffUsiem  model 
prediots  an  explosive  bubble  growth  at  only  .7®C  superheat*  and  affords 
no  asymptotic  solution  at  all  above  this* 


_ 

“ 


IV.  THE  DXNAMiC  PROBLEM 


The  eq[uatl<m  of  motion  of  tho  vapw  bubble  \^11  is  given  by 
6Qa  IX(68)y 


RR  + 


^ ss 
2 ” 


p (T)  - p 


q 


/ liq 


(1) 


iidiere  R denotes  the  radius  of  the  bubble  vail  at  time  ty  R the  radial 

Velocity  and  E the  radial  aceeleratioa*  The  initial  conditions  for  the 

bubbles  considered  here  will  be  that  the  bubble  starts  ftrom  rest  ^^ith 

radius  R • 
o 

R(0)  - R^y  R(0)  = 0,  (2) 


In  the  case  of  the  growing  vapor  bubble  y it  will  farther  be  assumed  that 
the  bubble  is  initially  in  (unstable)  equilibrium y in  tMch  case  the  radius 
beoaues  determined  by  the  equation  of  motion*  For  the  oollapsing  bubble  y 
equilibrium  eonditions  will  not  be  assumed  y so  that  the  initial  radius 
remains  in  this  case  arbitrary* 

The  s\irface  tension  parameter  and  the  density  ^ of  the  liquid 
will  be  assuaed  constant  y and  equal  to  their  values  at  the  initial  liquid 
teoqterature  T • The  external  pressure  will  also  be  assumed  constant* 

V 

The  equilibrium  vapor  pressure  of  the  liquid  P__(T)  at  the  tejo^rature  T 

0<J 

of  the  babble  wall  cannot  be  assumed  constant  y however  y since  it  is  the 
pressure  difference  p - p appearing  in  eq.  (1 ) which  supplies  the 
driving  force  for  the  expanding  bubble  y and  this  has  been  assumed  to  be 
in  Initial  equllltoium  with  the  surface  tension*  For  the  case  of  the  ooUapsing 
bubble y the  temperature  at  the  bubble  wall  rises  sharply  near  the  point  of 
collapse  y and  the  corresponding  rise  of  vapor  pressure  within  the  bubble 
may  be  expected  to  influence  the  rate  of  collapse* 

The  dependence  of  the  equilibrium  vapor  pressure  on  teaperature 
oan  be  taken  l^om  equllilnrium  vapor  pressure  tables  y so  that  p (T) 
say  be  assuaed  to  be  a knoim  Amotion  of  the  tesq>@rature  at  the  bubble  wall* 

The  aquation  of  motion  t!ien  becomes  determinate  when  the  bubble  wall 
teiqporature  is  specified  is  terms  of  the  parameters  of  bubble  growth  or 
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eoUapse.  This  spscifioation  vlH  hsi^a  ba  assumed  to  be  given  the 
ui^rturbed  oonveotiOQ  solution)  eq«  ni(26}t 

T(R,t)  J*  *1=^“. 

J J (7) 

' X 

vhsre»  by  eq«  111(3))  the  teiiq)erature  at  r = » is 


(q(0)  = 0). 


The  equilibrium  vapor  density  y-'^a  vapor 

pressure)  a knoun  function  of  the  tes^erature  of  the  babble  vail*  The 
Ibnotl^  q(t)  in  eq«  (4))  vhlch  represents  the  aeeuoulative  effect  o^^  thezml 
radiation  absorbed  by  the  liquid)  may  be  taken  to  be  a linear  function  of 
tine.  Ita  effect  Is  to  initiate  tiis  growth  of  the  equilibrium  bubble  (by 
raising  the  vapor  pressure))  but  its  infLuenoe  is  extremely  transitory  and 
the  term  will  be  neglected  once  the  bubble  growth  has  begun.  It  will  be 
omitted  from  the  equations  of  motion  for  the  collapsing  bubble.  The 
parameters  L)  k and  D («  k//  0^)  appearing  in  (3)  or  (4)  will  be  taken 
as  constant)  and  equtd.  to  their  values  at  tlie  initial  liquid  tenqperature  T^. 

The  error  inourred  by  the  neglect  of  the  variation  of  L)  9 
0)  eto.)  with  temperature  is  not  significant  in  the  ease  of  the  es^anding 
bubble)  beoatise  of  the  small  tesqjerature  variation  which  occurs  (essentially 
the  initial  superheat  T.  - T,  ).  The  error  involved  may  be  larger  for  the 
collapsing  bubble)  depending  on  the  initial  tes^erature  aM  the  initial 
bubble  radius ) but  is  not  as  serious  in  this  case  as  the  failure  of  the 
basic  assumptions  that  the  liquid  is  incompressible  and  that  the  babble 
Remains  spherical*  The  trend  of  the  physical  quantities  which  describe  the 
collapse  of  the  vapor  bubble  is  the  same)  whether  L)  D)  etc.)  msy  anr  not) 
however)  and  is  given  ooxxeotly  by  the  analysis  to  follow* 

In  terms  of  the  constants 


h;  - 

J Vr 


m 
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(iiMoh  have  ths  dimensions  800“^ » respeotl’Tely)*  ve  najr  define  a set  of 
djjaensicnless  Tariables 

o 


R (T) 

*'“2rll’»-Paqt*>’>  5°  J^iry 

/ ©q  o 

in  terms  of  '«d3ioh  the  system  of  equations  to  be  solved  becomes 

*-  2 


f 

> 


(6) 


(7) 


2,.^  p 

Jo  v/u  “ ▼ 


(8) 


(9) 


^ = 5 (T). 


The  initial  oonditions  for  (?)  are 


(10) 


0 * 1 * ^ n s=  0. 

The  physical  quantities  ve  eventually  vish  to  find  are  then  given  by 

t.i  r-^ , 

« Jo  ' 

Il(t)  = R 
o 


s> 

.r  r iiMfr 


and 


0 k 


>/rr 


(11) 


(12) 
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The  Eicpanding  Vapcr  Babble  > 

The  pure  'vapar  bubble  which  is  to  grow  in  a superheated  liquid  is 
assumed  at  some  stage  of  the  superheat  to  be  in  tuistable  equilibrium  usder 
the  effeots  of  surfaoe  tension » vapor  presage  and  external  pressure«  The 
bubble  growth  begins  as  a result  of  further  superheating  t ^ch  increases 
the  vapor  pressure  and  upsets  the  equilibrium.  The  condition  for  equilibrium 
at  the  time  of  release  of  the  bubble  t = 0 is  that  ^(0)  ~ R(0)  = 0*  and 
henoe  by  eq.  (1 ) that 

^ 0 

£q.  (13)  fixes  the  initial  radius  R of  the  bubble.  As  has  been  noted 

previously^  the  nucleus  £t<m  which  an  actual  bubble  grows  is  not  necessarily 

2 

spherical $ and  its  surface  energy  may  be  appreciably  less  than  4n<rR^} 
however } the  nucleus  from  which  an  actual  bubble  grows  and  the  free  spherical 
vapor  biibble  of  radius  R^  are  both  in  imstable  equilibrium  with  respect 
to  growth  at  the  temperature  T^  and  external  pressure  Table  II  (p.  6%) 

gives  a set  of  values  of  R^  for  various  superheat  temperatures  in  water 
at  an  external  pressure  of  1 atm.  From  the  definition  of  (^(T)f  eq.  (6)i 
or  the  differential  equation  (7),  the  equilibrium  condition  (13)  is 
equivalent  to  the  condition 

)if(T^)=-1.  (14) 


As  the  bubble  grows » the  tenperature  at  the  bubble  wall  decreases 
toward  the  boiling  point*  Inasmuch  as  liquids  will  ordinarily  support  only 
a few  degrees  of  superheat*  the  temperature  variation  involved  in  the  growth 
is  small*  and  an  approximate  expression  for  the  dependence  of  vapor  pressure 
on  temperatiare  will  suffice.  For  1 atm.*  a close  fit  to  equilibrium 

vapor  pressure  data  for  water  between  100®C  and  110®C  can,  be  obtained  ty 
taking 


p (T)  ~ p 

.it  ^ " =A(T-T^), 


(15) 


with  Tjjj  - 1C0®C  for  water*  and  A = 40*800  c.g.s.  unj.t8.  Ey  comMaiiyt  (15) 
vi*^  (S)s  (6)*  (9)  and  (13)  or  (14)*  me  obtains  for  fi  the  relation 


raaraQisenFRjcaKusniitdVj 


- 3 / — 


2 2 
% « ^ 


f 

R?  a^  J 


dv  M 


a)  dv 


The  term  involving  q in  eq.  (16)  is  extremely  snail,  and 
therefore  of  importance  only  for  a minute  portion  of  bubble  growthi  it  upsets 
the  initial  equilibrium,  For  a temperature  rise  of  1°C/nin  in  the  liquid, 
this  term  is  of  order  10  and  it  will  be  neglected  once  the  bubble  growth 
has  begun,  T*)  fix  the  model,  ve  shall  take 

~ q(t)  = at,  (17) 

correspvinding  (see  eq.  (4))  to  ten^erature  rise  of  1°C  in  see,  in  the 
liquid  far  from  the  bubble.  Then  from  eq,  (12), 


r -rf 

RT  k Jo 


where  the  constant 


1 = 


2 3* 
o 


In  kseping  with  the  above  (18)  may  be  approximated  by 


" q(t)  = "Xu, 

R~  a k 
o 


Because  of  the  small  tenperature  range  occurring  for  the  expanding 
bubble,  we  shall  further  approxi^t©  6 in  eq,  (16)  by  unity,*  and  write 
9q«  (16)  as 


- m = 1 + Ttu  -A  /“ 

^ Jo  V 


The  error  involved  her©  In  setting  ^ = 1 ( (T)  ^ ^ 

may  be  eatiniated  from  the  results  given  below  for  the  tesmerature  variatlt^. 
It  is  found  that  the  ratio  ^ remains  less  than  5 per  cent  at  any 
time  for  the  growing  bubbles' considered  here.  This  ratio  is  identically 
the  ratio  Q»f  vapor  velocity  to  Uquid  velocity  at  the  bubble  wall,  which 
has  been  discussed  previously. 


St 

I . 
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dz 


where  8’  (v)  aad 


At 


rJ  o^‘ 
0 


/ \ 


The  tsjretem  of  equatioua  for  the  growing  bubble  thus  simplifies  to 


l|U[J/3„2]  =1  -X„-^  J 


B»  S 0, 


- 1 at  u - 0, 


(23) 


A solution  to  eq*  (23)  will  be  found  in  four  papts*  eorrespcaiding 
to  four  (overlapping)  phases  of  bubble  growth*  which  may  be  labeled  vhs 
"relaxation  period" » "early  phase"*  "Intermediate  phase"  and  "asyn^stotic 
phase" • 


Since  the  bubble  growth  starts  from  equilibrium*  we  shall  put 


w 

K A * 


(24) 


and  assume  that  initla].ly  w(u)  and  its  derivatives  are  small*  On  neglecting 
the  second  or  higher  powers  of  w*  w'  * •••  » or  products  of  such  terms*  one 
may  rewrite  eq*  (23)  in  an  approximate  (linearized)  form  as 


W”(u)  - »(u)  = 3 lu  - 3/^f”  y'  (t)  4v  , 


/ t/o  — V 


w(0)  = w«  (0)  = 0* 

% putting  y(8)  = X.  taking  the  laplace  transform  of 

(25)  with  respect  to  u*  one  obtains  for  y(s)  the  equation 


(25) 


y(8)  - y(s)  =5  ^ - 3/x,sy(s)  • 
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„/g)  » ijk.  I ' 1 

8 8 - 1 + 3^/1V8 


(26) 


Id.  order  to  matoh  a later  solution,  we  shall  be  mainly  interested 
in  the  asymptotio  fora  ot*  the  8<^lution  (25)  as  u , obtainable 
CrosB  the  expansion  of  (26)  about  t}»  singularities  of  y{s)»  It  is  possible 
to  obtain  a solution  to  (25)  in  olicsed  form  by  the  means  described  below, 
and  also  to  write  a series  expansion  of  w(u)  in  powers  of  u firom  the 
lAurent  expansion  of  (26)  about  s = 0,  although  these  iiri.ll  not  be  used  here* 
The  roots  V's  , say,  of 

s^  + 3/a^  l/ffs  -1=0  (27) 

correspond  to  simple  poles  of  y(s)*  Eq*  (26)  may  therefore  be  expanded 
in  partial  fractinna  using  the  factors  iixdicated  in  (2‘/)*  For  a given  root 
\/s  = V^,  one  obtains  terms  of  the  fom 

_2 


± 


*^(  \/s  - v'^ ) 


_L 


v/s  - ^ 


(28) 


multiplied  by  constant  complex  coefficients*  By  the  use  of  the  Laplace 
inversion  integral  it  may  be  shown  that 


-1 


± 


J 

for  all  complex  \/^ , and  hence  that  for 


— ...X-..  + [“I  + erf('/j0u)]  (29) 


I < jarg  V^j  < rr, 

(29)  vanishes  as  u -*  « , It  follows  that  the  behavior  cf  w(u)  as  u » 
is  determined  by  thrjse  singularities  of  y(s)  for  which 


|arg  \/s|  < I* 


(30) 


Actually,  there  is  but  one  root  \/3  = of  iZl)  satisfying  condition  (30) 

for  0 < w , and  this  root  la  real. 
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Too  residue  of  3r(s)  at  a is  given  ty 


3X  -I 


2/^ 


or  since  ^ satisfies  (27 )» 


Hence  as  a , 


SO;  that  as  u -*  <»  ; 


3 X ^.S 
yS^  3^^  + 1 


r(s)  ~ 


w(u)  - 


^6  i3fi> 


£L 


+ 1) 


(31) 


^ (3^-^  + 1) 

where;  again; ^ in  (31)  is  that  root  of  eq»  (27)  satisf^rlng  condition  (30)# 
Alternatively;  eq«  (31)  xnay  be  written 


u ~ -s-  In 


p\^p  1/ 


* w 


w ■»  CO  , (32) 


^ince  the  transfcsm 
[T  - Tj  • 


r 


dv 


i_v/  o yu  - V 


\^s  y(s) 


is  asyiqptotic  to  y(s)  S'*-  ; it  follows  that 

^ ” ^00  ~ tt  -♦  00  , (33) 


Moreever;  to  the  degree  of  apprcximatlon  used  in  the  linearisation; 


u = at; 


w=^3(|--  D# 

o 


1 

J 


(34) 


II 


l^om  (31 ) f (33 ) thus  obtain  the  relations 

(35) 

(36) 

(37) 

(38) 

Thus  the  bubble  radius  remains  practically  equal  to  until  the  time 

tat  - vhen  it  begins  to  increase,  reaching  2R  at  about  t = t . 

0 o 

A tabulation  of  the  significant  parameters  in  e^s*  (36),  (38)  for  water 
at  1 atm.  will  be  found  in  Table  I,  with  the  choice  of  a = ,01®c/sec.* 


TABLCi  I. 

Parameters  of  the  Relaxation  Period 


R cm 

see 

ft 

l/a^  sec 

— -I  1, 1 

.01  t^ 

0 

102®C 

1,56  X 10”^ 

7.34  X lo"^ 

5.05  X 10*^ 

1.97 

7.34  X 10*^ 

104^G 

.75  X 10“^ 

8.08  X 10*^ 

. -6 

4.48  X 10^ 

3.30 

8.08  X lO”*’' 

106®C 

.48  X 10“^ 

3.09  X 10**^ 

1.56  X 10"^ 

3.71 

3.09  X ICT^ 

^8  choice  for  the  parameter  a correspo^s  roughly  to  the  rate  of 
the  temperature  rise  observed  by  Dergarabedlan(15}  in  his  experiments 
on  bubble  growth  in  superheated  water. 


R ~R, 


1 + 


2'X 


^.2 


+ 1) 


T ~ T^  + at  - 3"5  (r*** 


^c/t 


^Sr  defining  a time  t^  by 


one  may  write  eq.  (35)  as 


sj 

+ 1) 


O j(g(t“t  ) 

R - P.  [1  + e ' ® 1 . 


- 62  - 


Because  l/o^  « t^g  the  hvibble  growth  appears  to  start  abruptly  near 
rather  than  at  the  tline  o.f  release « t = 0* 

For  given  initial  eonditicms  T^g  p » the  duration 

C*  00 


t 


0 


= -t. 


In 


t 1) 

“TtC 


(39) 


of  the  relaxation  period  is  eoo^letely  determined  by  the  heat  source  function 
q(t)g  iae*  by  the  constant  ')L  % upon  which  it  depends  legarlthmioally« 
However j it  la  evident  frm  Table  I that  the  heat  scairee  term  (at)  in 
eq.  (36)  becomes  negligible  in  comparison  with  the  other  terns  near  the  end 
of  the  relaxation  period*  Brom  a physical  standpoint*  this  means  that  at 
later  tines  the  bubble  behavior  is  independent  of  the  rate  of  Increase  of 
stqperheat  vhieh  initiated  the  growth* 

The  asymptotic  formulas  for  the  linearized  eq*  (25)  presented  above 
are  accurate  over  a range  roughly  defined  by 


u 


p 


w « 1, 


or 


< t < t^ 
o 


Because  of  the  smallxiess  of  ^ * v(u)  increases  over  the  range  by  a factor 
of  several  powers  of  10  (about  10^)* 


written 

1 J- 

6 ® dw 


In  terms  of  v = In  z as  independent  variable*  eq.  (23)  may  be 
13w 


e ^ 


w 


= 1-0  ^ '’/x- 


j; 


eTdv 


'^0  \/u(w)  - u(v) 


(40) 


vrlth  the  neglect  of  the  heat  source  term*  For  small  w*  this  reduces  to 


} -i  f" 3s . 

6dvW  ■ 


/ 


vhleh  is  patisfled  by 

u » -J—  in(Kv)  C41 ) 

V 

for  arbitrasy  Ky  provided  that 


% 

Since  tl^  intep'al  has  the 
from  the  disoiisaloa  of  the 
various  roots  of  (42)  i the 
eondlticm  (30),  v = • 
ve  mast  further  set 


value  \/tt,  eq»  (42)  Is  identical  vlth  eq*  (27)* 
relaxation  period,  it  la  clear  “Ujat  of  the 
one  to  be  chosen  Is  that  one  vhioh  satisfies 
Xn  order  to  natch  the  previcus  solution,  eq«  (32), 


/?(3/^  + 1) 

~Tir“ 


(43) 


in  (4I  )•  ^ 

It  is  annarent  (41 ) that  the  derivative  ^ of  the  solution 
tt(w)  of  (40)  has  a simple  pole  at  w = 0,  vhioh  suggests  a solution  of 
the  form 


or 


du  ^ 1 
dv  ^ 


i 1 


jW  T SgW 


0, 


^2  2 

+ w +•  w + 


/- 


2/S 


1 1 1 \ 
\****/ 


(45) 


ly  abbstituting  (44)  into  the  integral  of  (40),  one  obtains* 


r - ; f -)  1 + I 12’/*  + 2(1  - 2^^)  a,] 

V o \A»(v)  “ b(v)  ' (. 

+ ^ [3^'^2+  6(2^/2_  31/2^^^  ^ 3(1-2*2^/2+  33/2)a2  ^ . 

(46) 


^ eq*  (40),  this  must  equal 

1 - a-*^  - J = a (1  .^2)  - ^ +^2(11  - «.,)) 

^ 2?  2('j-j§  ..  129a^  + 54a^  ^^2^^  4^  **»  » (47) 

* ' 1^'e  i^pe'odix  C for  the  avaluatlcxi  of  tiie  integral* 


The  paraneters  •••  are  ttisn  found  equating  the  ooeffioients 

of  corpeapfmding  powers  of  w in  (46),  (47)*  A tabulation  of  the  first 
seven  of  these*  for  varying  superheats*  will  be  found  in  Table  II  tot  the 
ease  of  water  at  1 atm*  external  pressure* 

The  time  oorresponding  to  (45)  becomes 

t as  1 r u'  (v)  dv 

® J o 


a -I. 


"^(3y6^  + 1) 

~TZ 


n 


w 


+ (a^-  |)  w 


^*2  "■  3 *1 


t|) 


2 


» 


(4S) 


the  logarithmic  tern  having  been  chosen  to  natch  eqs.  (32)  > (34) » The 
teiiqjerature  may  be  found  £tm  eq.  (46)* 


Asymptotic  Phase. 

During  the  early  phase  of  bubble  growth*  eharaotarlsed  the 
relaxation  of  surface  tension  * there  is  a rapid  rise  in  the  radial  velocity 
R of  the  bubble  wall  until  the  cooling  effect  of  evaporation  becoses 
important*  Ihe  rate  of  bubble  growth  thereafter  is  ecaitrolled  by  a bsi»nc-s 
between  the  rate  of  evaporation  and  the  rate  of  cooling  it  produces*''*’  Thus* 
wMle  the  vapor  oavily  grows  1^  evaporation  (since  the  vapor  velool'ly  Is 
ne^igible)*  the  motion  of  the  liquid  is  caused  by  the  difference  between 
internal  and  external  pressure*  Ho^^^ever*  an  increase  in  the  evaporation 
I'ate  tends  to  deci'ease  the  pressure  difference* 

•»'  " ' ' ' ' 

'Bie  effects  of  liquid  Ineirtia  are  Izoportant  in  de'bernining  the 
bubble  growth  near  the  tine  of  naxlraga  radial  accsleration*  In  the 
a^pmptotlo  phase  of  bubble  growth*  however*  the  inertial  term 


1 A, 

6 dis 


in  the  differential  equation  is  of  smaller  order  than  the  swfhee  tension 
1 This  point  will  bo  returned  to  below* 

^TTf  • 


TABLE  II 


Parameters  of  the  Early  Phase  of  Vapor  Bubble  CSrowth 


^0 

102 

103 

■*04 

1<^ 

106 

R X 10^  cm 
o 

1.558 

1.019 

.7507 

.5901 

.4832 

o X 10?  seo~^ 

1.797 

3.391 

5.356 

7.677 

1.035 

t ”0 

1.120 

1.023 

•9628 

.9205 

.8880 

.5598 

.3411 

.2407 

.1841 

.1480 

/ 

i .1101 

.2632 

.4168 

.5340 

.6177 

^1 

2.0915 

2.0322 

1.9763 

1.9456 

1 .9292 

^2 

2.1577 

2.0547 

1 .9463 

1.8852 

1.8526 

^3 

1.4633 

1.3761 

1 .2732 

1.2129 

1.1807 

.7359 

.6864 

.6224 

.5831 

.5620 

A 

”5 

.2946 

,2722 

.2427 

.2236 

,2132 

.09858 

.08945 

.07886 

.07143 

.06726 

.02857 

.02510 

.02205 

.01960 

.01819 

66 


¥ 


Notff  it  is  elear  that  th@  bubble  snxf9t  contiime  to  slues  a 

stationary  bubble  Is  at  the  temperature  of  liquid  superheat  and  therefore 
has  a high  internal  pressure « Hgnee»  the  ttanperature  at  the  bubble  vail 
ssust  continue  to  drop  betsause  of  o'^nperatieais  But  the  tenq>erature  cannot 
dr<^  h&lmn  the  boiling  point  and  still  maintain  the  pressure  difference 
necessary  for  growth*  It  fdllnus  that  the  teoperature  of  the  bubble  wall 
xmiat  approach  a limit  as  t ; and  this  fact  is  sufficient  to  characterise 
the  asysptotic  phase  of  bubble  growth*^ 

It  is  perhaps  worthwhile  to  demonstrate  at  this  point  that  the 
limiting  temperature  predicted  by  the  mathematical  model  is  what  one  would 
e3q>eot  m physical  grounds  - the  boiling  teiiperatare  of  the  liquid  at 
the  external  pressure » for  the  sake  of  conslstanoy  and  to  Justil^  statements 
made  in  section  III  abo^e*  The  differential  equation  i with  the  neglect  of 
the  heat  source  term*  may  be  %fritten 

If  tile  last  term  on  the  right  (the  inertial  term**)  tends  to  vanish  as 
u *»  « f s 00  * then  the  limiting  value  of  the  temperattar©  integral  on  the 
left  side  of  eq*  (49)  is  unity: 


/^X  "*1  as  u -►  «> , 

vu  - y 

^ actual  teoperature  is  given  by 


T=T„-f 

V A 01  - V 


30  that  by  (50)* 


o vu  - V 


T^  - T ~ -r“  as  u • 

O 


The  slew  teoperature  rise  due  to  irradiation  is  here  neglected< 
The  kinetic  energy  of  the  liquid  is  given  by 

f (J  dr  = a«^. 

Ojj  *i/  yJ  O 


(50) 


(51) 
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According  to  'Uie  definitions  (22)  arid  (5/j 

i 


^ - 2<r- 

" A * 


and  80  frtm  the  definition  (13)  of  R^» 

yu->  A y/£> 

Ey  eq*  (15);  this  states  that 

i 

S T - T. 

^ 0 b 

(y^  being  considered  a constant) i and  therefore!  by  ooo^arlson  vith  eq®  (51 }» 

T " T.  as  u oo  • (52) 

b 

The  conclusion  (52)  depends  on  the  equation  of  motion  only  to  the 
that  It  follows  S^m  the  aaya^totic  vanishing  of  the  inertial  term. 
One  can  easily  shsvt  however,  that  tl®  inertial  term  mat  vanish  tdiether 
(52)  holds  or  not,  provided  only  that  the  temperature  approaches  some 
limit  as  u •♦  «>  • For  this  implies  that 


8^  dv 

"■  ' ~ const, 

o /u  - V 


as 


CO 


(53) 


Idiltlplieation  of  (53)  by  and  integration  from  u s o to  x 

V*  - u 

yields,  alter  a change  in  the  order  of  integration, 

X 


r*  = [’'..wdT  f 

VO  U J'H  /u  - V ^ O 


V s/^ - U l/u  - V 


= ff[s(x)  - 11,  const.  • 2,<i^  (54) 

The  vanishing  of  the  inertial  term  in  (49)  then  follows  ITom  eq.  (54),  idilch 

shows  in  feet  that  the  inertial  term  is  of  a smaller  order  of  magnitude  in 

-‘1/3  — 

the  aqyggjtotio  range  than  the  svirface  tension  term  s ■ The  constant 
in  (54)  is  1^  according  to  eq.  (50),  so  that  eq.  (54)  yields 

* (tt)  ” u *•  <»  « (55 ) 

!f^4^ 
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Eq#  (55)  deoeribes  tl»  asnmtotle  bubble  gr<^hf  but-  Is  not  yet 
taisef^iLfy  sinoe  it  provides  no  means  of  znatehlng  the  indicated  ssysq}totlo 
solution  cf  eq«  (49)  vlth  earlier  sclutlcms.  The  possibility  of  matehing 
solutions  depends  on  the  possibility  of  shifting  the  asyn^totio  solution 
In  t (er  in  u)  so  as  to  account  for  the  relaxati^  period  of  bubble 
growth*  It  is  necessary  that  one  be  free  to  shift  the  asynptotio  solution 
since  the  duration  of  the  relaxation  period  depends  completely  on  the  choice 
of  tl%  heat  source  termy  while  the  subsequent  behavj.or  of  the  bubble  is 


independent  of  this  term* 

The  means  for  makir.\g  an  arbitra:^  time  shift  is  ftirnished  by  noting 
that*  in  addition  to  the  asys^totic  form  of  solution  (55 )»  eq*  (49)  also 
possesses  the  solution  s(u)  » 1*  Tht;^  the  complete  asymptotic  solution 
may  be  described  by 

b(\i)  =s  1 y 0 < U < Uj  I 1 


> (56) 


= 1 

Vo  )/xi  - T 
&qm  eq*  (12)  y the  time 


jL  [^7/3 

A*  I 


1/3  6a*  du 


corresponding  to  eqs* 


a*^]  y u > u^ . 
(56)  becomes 


t 


/: 


175 


§SL 


(v) 


so  that  u^/a  here  represents  the  diaratlon  of  the  relaxati«m  period*  The 
time  shift  may  be  introduced  explicitly  into  the  asymptotic  solution  by 
using  the  fbot  that  if  a(u)  is  a solution  of  eq*  (56) » a(u  t u^)  is 
alse  a solutiony  with  delay  period  (u^  ~ 

A consistent  scheme  for  continuing  the  asymptotic  solution  may  be 
found  by  takiug  ths  solution  to  be  of  the  form 


a(u)  = 1, 

a(u)  "*  f I + ^ 


0 < u < u*  # 


i 


b^  la(uHa  ) / ^ ® 


> (58) 


to  floven  toms*  ^ere  is  a constant**  Vhen  the  coefficients  b,,  haire 
been  determined*  the  difference  (u^  - vl^)  is  fixed  by  the  requirement 
B (o^ ) - 1 • The  matching  vith  earlier  solutions  is  then  ecconiplished  by 
adjU8'UJ3g 

When  (5d)  is  substituted  into  the  integral  of  (56)  the  result  is'*’ 


r sjLk 

Uj  s/n  - 


1 


4-  .w^nA, 


1 + sS9266  ^ < 5 


(u-u^j 


2/6 


b. 


(tt-U^) 


3/6 


+ .47545 


(u-u^) 


475 


+ .27450  2 

(i4-u 

c 


u - u. 


(59) 


(^) 


By  eq.  (56)*  the  eaqiressian  (59)  is  also  a^snptotie  to 

This  may  be  expanded  by  (58)  to  give,  on  equating  coefficients  of  correspond- 
ing poysrs  of  In  (59),  (60),  a set  of  successive  equations  for  the 

parameters  b^ , . « » , At  each  step  one  has  (as  for  the  early  phase 
coefficients)  a linear  equation  for  the  iinknown  parameter.  A tabulation 
of  these  parameters  for  vax^ious  superheat  teu^eratures  in  water  at  an 
external  pressure  of  1 atm.  is  given  in  Table  III. 

leading  terms  in  the  asymptotic  solution  are 

b = ~ ^ f ■ * J • 


O 


o 


1 + 0(u”^/'^)  J 


Hghor  terms  are  of  the  form  [in  (u-«  )]**/(u*^a . where  n and 
k are  integers.  ^ 


The  natch  is  better  obtained  In  praotlco  by  shlfi;ing  the  asynqptotlo 
R(t)  ourvsa. 


Sea  Appendix  D for  the  evaluation  of  the  integral 
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TABL2  III 


Parameters  of  the  A^nnptotlo  Phase  of  Vapor  Babble  Chrovith* 


T ®C 

0 

102 

103 

104 

105 

106 

H X 1(r  cm 
o 

1.558 

1.019 

.7507 

.5901 

.4832 

a X 10*^  sec”"^ 

1o797 

3.391 

5.356 

7.677 

1 .035 

t °c 

1.120 

1.023 

.9628 

.c«205 

•8880 

.5593 

^3411 

.2407 

.1841 

.1480 

‘>1 

“1.073 

- .9099 

- .8101 

- .7409 

- .6890 

^2 

- .4709 

- .4624 

- .7122 

-1.334 

-2.506 

-3  ^ 

- .2339 

-1.481 

- .3586 

- .9725 

-1 .972 

» .5534 

1.258 

- .1970 

"4.598 

-2.510 

”3 

"1.775 

3.064 

4.298 

18.50 

65 .84 

\ 

" .7423 

- .2670 

.6166 

- .4736 

-48.65 

1 
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Tboa 


+ o(r''^)J  1 


n 


2 t.2 

a R~ 


^ ^ Ms 


fir  t&  ttamo  of  the  erlglnfiLl  pkQraloal  eonatsntSf 


(61) 


-*yi 


’‘K  - \) 


v^r<. 


i 1 4=  0(t  "*) J 


T - (T^j  - ^ ^ 0(t“^/^)J  , 


(62) 


Ths  es^^Biptotio  ralatlcn  has  been  dleeusaed  pre^^ly*  The 

radius  jralatioo  has  also  a sinqple  pbgrsical  interpretation » which  may  be 
h^»  Ey  differentiatir^  ^ first  of  equations  (62)9  one  et^tains 

(T.-Tj  j 

itrlT  ’ 


S -/| . 


k 


s/Bt 


y. 


eq'  O' 


vM^  i$ay  also  be  written 


A»B®  • k'  ^ 


" -% 


In  ^s  foray  one  oay  readily  recognise  the  heat  transfer  relation  holding 
at  the  bubble  vail  f the  ri^t  side  gi'ving  the  heat  gain  in  the  -vapor  and 
the,  left  side  giving  the  heat  loss  zron  the  liquid*  The  teiqperature  gradient 
is  here  es^ressed  in  terse  of  ton  ratio  of  the  tenpsrat'jire  drop  oeeurring 
aaysqptotioally  at  the  bubble  ^11  to  tl^  thickness  of  the  thermal  boundary 


layer  in  iMoh  it  occurs*  Ibe  particular  choice  / 3 ^ the  eharaoterisiio 
diffusion  length  could  not»  of  course*  have  be^  predicted  beforehand* 

Uiile  the  leading  terns  of  the  asjraqjiotio  expansions  ^ven  in 
(62)  shbV  the  essential  -variation  of  the  physical  quantities  ^ch  describe 
the  bubble  grcw^t  they  are  of  liaited  usexUlneas»  and  nay  be  in  error 
as  sueh  as  40  per  oent  (depending  on  the  superheat)  idiiXe  ths  bubble  radius 
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Is  still  aaalls  accurate  esqirdssious  may  be  found  by  carrying  out 

the  tine  intsgratlw  indicated  in  eq*  (57)  and  subiitituting  in  the 
coefficients  f^om  Sable  III*  The  result  of  such  calculations  is  presented 


in  lig»  3 1 whioh  follows  asysiptotia  sdutions  down  to  a radius  of 
4 X loT^  m.  for  the  indicated  sttperheats  in  water  at  1 ata*  external  presstire- 
Inasmuch  as  the  duration  of  the  relaxation  period  ntay  oe  chosen  arbitrarily 
so  far  as  these  a^p^totlo  sc^utions  are  oonoemed*  the  tins  scale  is 
detersdned  only  to  within  an  arbitrary  additive  constant  (the  constant 
u^/a  mentioned  above)  whioh  stay  vary  ftca  one  curve  to  another*  The 
actual  spacing  of  the  curves  as  presented  was  chosen  sa  that  the  tine  inters 
oepts  at  R ~ .004  cm  were  equally  spaced* 

The  experimental  evidence  available  thus  far  covers  only  the 
asymptotio  range  of  bubble  growth*  Observations  on  the  growth  of  vapor 
bubbles  in  water  have  been  made  by  Dergawthedlan^^^^^  and  are  presented 
in  Elgs*  4f5»6»  together  with  the  theoretical  predictions*  The  tbsereWeal 
curves  were  obtained  hy  graphical  interpolation  from  the  set  of  eur'Tes 
plotted  in  Fig*  3*  The  tins  origins  both  the  the^etical  curves  and 
the  csqperimental  points  are  arbitrary!  so  that  a time  translation  of  the 
theoretical  curve  has  been  made  in  each  case  bo  give  the  best  fit*  The 
agreement  is  seeny  however*  to  be  very  good* 

The  lag)Ortanae  of  the  heat  transfer  at  the  bubble  wall  is  shown 
in  Tig*  4!  where  the  theoretioal  curve  obtained  with  the  neglect  of  this 
effect  is  also  plotted**  The  asyn^totio  form  of  tills  solution  is  readily 
obtained  firom  eq*  (1)  by  setting  (T)  = P«^(T^)  there.  The  differential 

^ ^ w • 0q  • gq  0 

e^iatioQ  say  be  written 


with  the  help  of  the  definition  (13) » and  yields  on  integration  frm 
to  Rf  to  Rf 


*2 


( “f  ) + 


JLO'' 


The  solution  for  the  motlan  of  a bubble  under  ooastant  vapor  pressure 
conditions  was  given  by  Rayleie^('^)  and  applied  by  him  to  the  case  of  a 
oodlaiMlng  bubble* 


R X 10  CM 


Fig.  3 - Asymptotic  radius  versus  time  curves  calculated 
for  water  at  1 atm.  external  pressure  and  the 

Indicated  superheat  temperature.  : 

j 
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t- to-M!LLISEC. 

Fig.  5 - Comparison,  of  theoretical  radius-time  values  with 
three  sets  of  experimental  values  obtained  in  water 
at  1 atm.  external  pressure,  superheated  to  104.  5°C. 


As  R «x>  , this  gives 


m 


R 


$ 


which  is  a oonatant* 

Ebqperinents  have  recently  been  performed  Dargarabedian  on  vapor 
bubble  growth  in  pure  rates  of  bubble  growth  in  this  liquid  are 

o«!55arsd  with  those  of  water  ut  the  same  value  of  the  temperature  difference 
(T^  •"  T|j)  and  at  oon^jarable  tlEioSf  they  should  bo  about  In  the  aano  ratio  as 
k/L for  the  two  liquids.  This  constant  is  3*5  times  greater  in 
water  than  in  carbon  tetraohleride.  Dergarabedisn's  bbaervationa  (ua  bubble 
growths  are  in  good  agreement  with  this  value. 


Intergadiatw  Solution^ 

The  early  phase  and  asymptotic  solutions  presented  above  join 
in  the  neighborhood  of  the  wfl.'ytnium  valtte  of  dz/du«  While  the  match  of 
these  solutions  is  fairly  good  for  all  superheats » it  is  nevertheless 
desirable  to  have  available  a solutioa  which  eovera  the  critioal  region » 
to  facilitate  the  mtohing  process. 

The  intermediate  solution  presented  here  idll  be  an  esq^ansion 
about  the  point  u » u^  defined  by 

Since  the  early  phase  solution  is  not  asstmied  to  be  accurate  at  this  point f 

the  actual  value  of  u^^  or  a(Uj^)  is  unknowno  In  order  to  determine  these 

quantS.tieSf  we  require  that  the  intemsdiato  solution  and  its  derivative 

shall  coincide  with  values  obtained  Arcm  the  early  phase  solution  at  a 

point  u ~ u idiere)  that  solution  is  accurate.  The  eiqpansion  about  the 
e 

inflection  point  u^  is  constructed  as  follows 

The  dif Acuities  ^diiob  arise  for  any  such  intermediate  solution  are 
connected  with  the  exi'ansion  of  the  convolution  integral  in  the  differential 
equaid-on.  Thas»  a solution  about  a knoi^  point  (such  as  Ug)  which 
assumes  the  integral  and  all  of  its  derivatives  to  be  known  does  not 
aotually  make  use  of  the  information  given  ly  the  differential  equation » 
while  a solution  about  a known  point  which  uses  the  minimum  of  data  necessary 
(the  valiie  of  s assd  say,  at  the  point  of  expansion)  involves  about 
the  same  procedure  as  that  given  below. 


k&Bwm  that  s(u)  has  an  expansicsn  about 

8(tt^  + x)  * Z[1  + c^x  + -^  + •••  1 (2  « z(u^)), 

so  that  2 

z*(u^  + x)  ~ Z[o^  + ^ + •••  1» 

z"(Uj^  + x)  = Z[c^x  + •••]  • 


(64) 


The  expanaicua  (64)  are  to  be  substituted  into  th€^  differential  equation 
(vith  the  heat  source  term  omitted) 


i 1-1)5 -^r 

[_  a*'*'  ' </©  /u  - V 

HoW}  the  integral  in  (65)$  evaluated  at  u + x»  ia 

'U+X 


pi*  aMz2^.,.= 

p — zLCsI-to  + f 

)o  \/tt  4 X - V 

*^0  /U  + X ~ V J 

(65) 


s/n 


<Jo  \/a  + X - r '^0 


a’  (u  + X - v) 
0 


= f + f £ <u  - v)  fs- 

O ‘/u  + X - V ks©  Jo  s/v 


(66) 


v&li.d  for  sufficiently  amall  x$  But  for  small  X(  z'(x)  s 0,  so  that 
when  (66)  is  valid,  the  first  integral  in  (66)  is  negli,gible*  Thus  for 


«U.4x 

^ g’tV 


X - V 


M T 4*  X X ^ T ^^***  4'  * * ' 

~ ^0  n ^2  zi 


as^,  where 


T - dr 

'1'  J«  vq — ‘ 


(67) 


(«8) 


- V 
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If  «.  were  knovu*  «9  could  ua©  (67)  ±n  eq*  (65)  to  obtain}  cm  equating 

«L 

eoeffiolenta  of  eorreaponding  powers  of  x in  (65 ) f a set  of  relatiais 


^ (69) 

®3  “ °1  •••  » J 

to  solve  for  the  parameters  c^j  o^*  •••  • However#  the  are  not  yet 
detex^ned* 


Since  the  solution  to  (65)  is  assumed  known  for  u < u^# 


let 

_ dv 

9 

(70) 

K 

W O '/uj^  “ V 

“k 

\/a£  - V 

9 

(71) 

so  that 

+ 

ii 

(72) 

in  (68)#  and  set 

J = - 

%•  c«u^-v 

(73) 

in  the  integralse 

Then  in 

L » for  instance# 
c 

n«  (v)  = 

a’  (u^  ” c)  - Z(o^  + ^ 

0^  + •••)» 

giving 

zlSsJU^ 

= 2Z  i/r  (e.  + ^ 

c)'^  + — ) . 

Thus 

^o"  ^0^*^  ) + 2Z  + 10  ®3'^^ 


- J^iS)  4-  2Z  /T  {-  I + “•)#  ... 


and  simll^ly 


> 


(74) 
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itay  ba  found  if  Z,  e^,  c^t  , and  are  known.  The  feasibility  of 

this  procedure  depenoa  cm  the  fact  that  a good  eatinate  of  o is  alrsa^ 
available  ft*catt  the  sar’iy  phaae  solution*  and  ttet  the 

yj)  = / * » 

'^O  + ^ - V 

are  slowly  varying  ftinctiona  of  ^ . This  follows  from  the  relaticsn 

« V<r>  = , (75) 

since  In  the  early  pl^ae  « a^*^\ 

Assuming  that  the  cssitributions  to  Srcsn  higher  powers 

of  S than  those  written  in  (?4)  s®y  be  neglected*  we  teradnate  the 
e:q>ansicaia  as  written  and  substitute  -Wiem  into  (69),  Together  with  the 
oonditions 

, 

~ *"  0^  d ^ Oj  o ^)» 

« ZCc^  +1 

firoD  (64)*  the  equations  y76) 

®3  ’ ^ f °i  ■ + azVS  (- 1 oj^ )] 

+ ^e^[J^(S)  + 2Z  (e^  -r  ^ OjS  J- 

constitute  a system  of  four  simultaneoUB  equations  far-  the  four  r.nt-n.-.^fna 
Z*  e^*  cy  S • Inasmuch  as  one  and  only  one  point  of  inflection  of  a(u) 
is  known  to  exist*  these  eqfoations  have  a xinique  solution. 

It  should  b©  noted  that  because  of  the  definitions  of  u and 
Bs  tJi©  waviimmi  radial  veloelty  #.  of  the  bubble  *^11  does  not  occur  at 
the  same  value  of  u (or  t)  as  the  maximum  value  of  s»(u)e  The  dlsoi'epaney 
Is  not  great  for  small  stqjerheats,  but  for  larger  mip9vha&%&  the  point  defined 
by  S ss  0 aomas  «^tQ  the  a^a^totio  imd  of  sCu)  curve. 


. -T 
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Soltitlcr*  for  the  Eicoatidlng  Bttbbla^. 

Die  results  of  the  above  theory  vater  at  Initial  teo^eratures 
-Q  ^ 106^0  and  external  pressure  1 atm*  have  been  plotted  in 

Sigs*  ?*  8t  and  9t  10*  The  bubble  was  taken  to  be  in  equilibrium  at 
time  t = 0 when  the  heat  source  q(t)  is  introduoeds  q(t)  was 
arbitrarily  ohosen  to  correspond  to  a ten^erature  rise  of  in  100  sec 
in  water  fhr  l^om  the  bubble** 


* 


See  the  footnote  on  page  61* 


“CM/SEC 


t - MILLISEC. 

Fig.  7 - Theoreiicai  radius  and  radial  velocity  curves  for  the 
growth  of  a pure  vapor  bubble  in  water  at  1 atm  ex- 
ternal pressure,  superheated  to  103°C. 
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Fig.  8 - Theoretical  radius  and  bubble  wall  temperature  curves 
for  the  103®  vapor  bubble  of  Fig.  7. 
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Figi  9 - Theoretical  radius  and  radial  velocity  curves  for 

the  growth  of  a pure  vapor  bubble  in  water  at  1 atm, 
external  pressure,  superheated  to  106  C. 
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Fig,  10  - Theoretical  radius  and  bubble  wall  temperature 
curves  for  the  106®  vapor  bubble  of  Fig.  9. 
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Tba  Gollapslne  Vapor  Babble* 

Baoauae  of  the  large  tes^r&ture  varlatiofna  tfhioh  occur  vfaen  a 
vepw  bubble  collapaes  in  vater  belcv  the  boiling  tes^raturej  a siiaple 
analytic  expression  for  the  vapor  pressure  or  vapw  density  variations 
cannot  be  found#  If  take  our  data  trm  equilibrium  vapor  pressure  and 
density  tables,  ve  conmit  the  treatment  of  the  problem  to  a numerical  one 
tffoa  the  beginning# 

^ system  of  equations  to  be  solved*  eqs*  (7)"(11)s  ia  unchanged# 
In  this  ease*  however*  the  vapor  pressure  at  the  Initial  tessperature 
is  less  than  the  external  pressure*  and  initial  oonditlons  of  dynamio 
equilibrium  cannot  prevail  for  the  p^u^  vapor  bubble#  There  is  therefore 
no  need  to  retain  the  heat  souree  term  in  eq.  (9)*  and  ve  shall  put 
q(t)  ~ Oe  We  continue  to  assume  that  initially  the  vapor  bubble  and 
surrounding  liquid  are  in  thermal  equilibrium  at  tea^rature  T^# 

It  is  eonvenient  to  transform  the  ten^erature  equation#  By 

# #i»^ 

multiplying  eq«  (9)  by  (x  « u)  *'"'  and  integrating  it  from  u = 0 to  x» 
one  obtains*  after  an  integratim  ^ parts,  the  relation 


^(u)  z(tt)  - 1 ~ ^ 0*(v)  v/u'-  v dv, 


where 


= T - 


The  system  of  equations  to  be  solved  beomee 

1 ^ ^ « 0, 

8 ' 

0 ^ 0{9)» 

^8-  1 =*-^  dv. 


(77) 


(78) 


8(0)  “ 1, 


9(0)  =0, 
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'Jhs  ^stea  (78)  vas  solved  immerloally  for  initial  tenqperature 
X s 22®C,  oxtazTial  pressure  p_=s  *544  atm,,  and  initial  bubble  radius 
(idiioh  is  undetermined  for  the  non-equilibrium  bubble)  = ,25  om.  The 
method  of  solution  is  given  In  %pendlx  S«  The  partieular  initial  data 
chosen  here  oorresp^^  to  values  idtieh  have  been  obtained  e3q>erimentallyk''*' 
Although  the  teoqperature  effects  become  significant  during  the 
collapse,  the  dynamics  of  the  particular  bubble  considered  l^re  differs  very 
little  from  that  predicted  by  the  Baylei^  solution  of  th®  problem  over 
most  of  the  collapse.  The  Rayleigh  solution,  which  neglects  heat  transfer 
effects,  is  readily  obtainable  £rm  (78)  under  the  assumption  that  0 
stays  constant,  and  equal  to  ^^(1^),  The  equation  of  motion  is 


Wiich  yields 


s(0)  * 1 , 2* (0)  = 0, 


^ . a)  + I ^ jjS/Sj 


(79) 


on  integration.  Since 

«2x  10^ 


is  much  larger  than  3/2,  eq.  (79)  may  be  appraxtmated  by 


or 

tdiere  the  negative  square  root  of  must  be  chosen  to  cnrrespond  to  the 
collapsing  bubble,  Ftok  eq.  (12), 

g-  = a z'  = - . ,1/'^  (1  - z)'/*, 

which  yields 

t=— i^r’ 

^-10  scxpciriasnts , perfonoed  at  the  It^dr  odTuamics  Laboratory  of  .the 
Institute  ef  Technology,  wore  reported  by  I'hS,  Flesset  in 
ref,  (13),  j 


d 


f 


A 


(hA„)- 


(80) 


y: 
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ou  i&t6gri!iti<»L*  The  time  oorresptmdi&g  to  the  lUll  ^tem  (78)  was  fouzid 
by  a numerical  integration  of  the  relation 


ufiing  the  values  of  z obtained  f^om  the  numerical  solution  of  (78)« 

A oon^arison  of  the  tvo  solutions  for  the  collapsing  bubble  is 

given  in  ilg*  11*  The  magnitude  of  the  radial  veloeity  of  the  bubble  wall 

obtained  lyon  the  numerical  solution  is  plotted  in  Fig*  12 » and  the 

oorresponding  temperature  at  the  bubble  vail  in  Fig*  13*  The  numerloal 

soluvicn  was  not  oarried  out  farther  than  shewn  in  Figs*  12  or  13 » because 

of  the  breakdown  of  the  assuaptims  underlying  the  theory  presented  here: 

the  parameters  L»  Dj  • etc*  begin  to  vary  significantly  near  the  end 

/ 

of  collapse » the  llqu5.d  veloeity  becomes  so  large  that  conq^resslbility 
effects  may  beooioe  important  in  the  liquid » and  the  spherical  bubble  shape 
boesnes  unstable  to  small  dlstoai^lng  influences* 


Figt  1 1 - CoTnparison  of  the  numerical  solution  of  the  text  (which 
” includes  heat  transfer  effects)  with  the  Rayleigh  solution 

(which  neglects  heat  transfer  effects)  for  a vapor  bubble 
of  initial  radius  . 25  cm,  collapsing  in  water  at  22®C  and 
an  external  pressure  . 544  atm. 
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The  taiblef  bel6w  give  repi^sent&tive  ^'aluas  taksii  E«N»  Dafsej^ 
Prop^x^ss  of  drdiiuu^  Ijatffs^tibs-Usi^a  PaibH^Msg  CorpOi:^tic^$ 

mOv  fork;  $ 1940)  • Vsluss  fss*  the  vapox*  are  tsheartain  C8xperljneatally)ia  th© 
third  Vlghifioant  ^guro  aiid  are  Ri3m\ih0.t  de^ifideat  on  pressUsrei  those  cited 
v^nrespos^d  to  pi.*@@syre8  of  1 ateie  or  beXov*  V&Xum  for  the  Ulqiild  have  been 
, !^feltr-^,13y  lifflitsd  to  fewr  aignifioact  ^garss<» 


T 


■'^kls^ior^’' 


• see 


«.  Viin 

9 eeo  • C 


}D;  (»  e^)  eaVa^c 


Water  Vapor- 
0 

•884 

1.39 


(ks  eccata^t  per  9ik»  o|*  vater  vaport 

B » n386,.*  1£^'« 


2.38 

1.41 

5*98 

•2^ 


e„ 

» C, 


1 #44 
2.68 
1 •^v2 


■7}  1 V -S^CSi''  ® ‘^0  : 

10*^;  sTs/m  * :Sfr  “®'’ 


fi-  X 1©  ^ • **C 

vt  /<.« 

f:s',1#  fsa^w 

■■  '■•A''  , s 

P ds^i/bia?^ 

*‘9  . 
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gyaluatltm  of  the  mffiiajnn  latssrals  (ifTa). 

In  ternci  of  the  nev  variable 


t = J-=-v£  , 
1 + 


the  integred  III  (4%) 
t 


la)  - 1 e 

Jo  JT^ 


iteL_; . «.  I-'/*  f 

J 


baeones 

1(A)  * 2 A 


>1!  t 


l=i«  .-aS 


-/: 


- -Ai 

J:±^^  t jfli 


© (1+t)^  ^ Jo  (I4t)^ 

l^'he  farther  substltutioin  t ^ 1/x  in  the  aeoond  inte^al  gives 

.1  . , -4! 

1 


'O  (I4t) 


,2 


e 


t JSl^  I 
/t  ‘/I 


i^Q> 


-J?lU 


(Ifx)' 


/s 


and  therefore 


2A 


I (A 


r 


o (1H) 


ldt_  e-A  Jt 

^ * 

V w 


= [1  + —=^3-]  / 

d(A^} 


(Itt)'  \/b 


Since 


X 


r«o 

Jo 


e 


-(1+t)x 


X dXf 


(m)‘‘  '-^o 

the  integral  in  (d)  any  also  be  written 


r.-A-^  r 

VO  /t  J 5 


'C& 


^ o 


/. 


(a) 


(b) 


(0) 


id) 


*“  CxfA^)t  ^ 


e 


v/t 


r*  y 

* I o**  X /--2- 

‘^o  V * + A 


2 


(e) 
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aftax’  a ohange  in  the  of  Integration*  Bgr  putting  y ^ one 

obtafjDS  after  bs&  integiration  by  parts 


V^*  f 


00 


6-  -—5^  = y/v  f"  (y->2)  *- 


- >r  r>a.  1 f 

- V”  tA-r  « e j 


_-y 

° 


■^2  “ 

>r 


= A v/^  + (1  ~ 2 X^)  o ^ ^ e"^  . 

x"  ^ 

!Che  us@  of  (f } for  the  intagr-al  in  (d)  then  gives 


(f) 


^«2X/?-2XV 


.>9 


or 


I(X)  “ 4X^  \A*  ~ 4w  erfo(X)» 

The  asya^totlo  formulas  f^  l(X)  foUov  f^oa  the  relations 

.2 


a 


^ ex'fcCX)  ~ ~::7^~7~  [t  - 


;i/n 


2A‘^  4X 


4 


h X 


CO 


arfo(X)  "•  1 - [1  - -.  ••»  3, 

yw  “ • • 


X-0* 


The  integral  in(47e> 


!(»=>  ['  --^ 

»'o  ™ 


X 

(. 


. 1+*^ 


(g) 


(h) 


am. 

/ » \/l  +/J 


J 


(i) 
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mj  bs  transfonaed  the  substitution 


/L^ 

1 + v/x 


(J) 


tu 


r1  1^2  .2+2 

lOO  = 4>  / ^2  e”  ^ ^ dt 


r [>t  + -^-^1 

* _t ^(1*^  ) 


p(  erfc 

Fop  ^ « 1 , t < /t  , 


At  + 


A(1  « t^) 


At  + 


]A(1  “ t^) 

and  the  brace  In  (k)  beconea 

v2 


= “+  o(A^), 


(k) 


1 “ \/ir  (^)  erfo(^)  + 0(A^)* 


(1) 


As  t in  (k)  increases  beyond  ^ ^ the  brace  in  (k)  drqps  rapidly  to  zero 
(A«  1).  But  £i*GB  eq,  (h)f 

1 - (^)  erfe(^)  ~2  ^ t ■*  1,  A 

80  that  for  small  A the  approximate  eaqireaaion  (1)  differs  frou  the 
braoa  la  (k)  only  in  0(  A ) the  full  range  of  t.  Since  the  factor 
outside  the  brace  is  unity  to  0(X'^)  vhen  t & A > apprcndmation 
to  the  integral  ^dsich  is  talid  to  terms  of  relative  order  A*'  is 

AX  dt  - V5  (i)  erfo 

= AA  - ZX®  ■T 9rfe(i)  ~ 1 + ^ 

(.  A Xv^  J 
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by  (h)  agal»« 


lO.)  *•  2 \/S  [1  “ 0(  X^)1  as  X 0. 

v/n 


(a) 


'^2 

fw  large  A the  brace  tn  (k)  remains  near  unity  until  t “ 1 =*  0(l/X  )» 
so  that  the  dominant  ihotor  in  the  integral  is  e**  ^ "t’  . It  Is  ocnvenientf 
therefor®  s to  put 


1 -t" 


(a) 


in  {k)f  and  %i»*lte  the  Integral  as 

r'  (1-t^)  .->iV 


i(»  = « / 


n 


0 (m^)^ 


\2.2 


2 1 
dt  1 “ ^ e*  erfo(x)  ^ , 


(o) 


In  the  region  ^re  e “ is  still  large*  a is  of  order  l/X 
and  X is  of  order  usdty.  Thus  in  tMa  region*  ve  may  expand 


ee 


erfo(x)  = erfo(  At  + ~)  « erfeCXt)  - Y 

^ \/v  ]s^ 


C(i) 


lC”1 

e J 


ki 


.2.2 


= erfo(Xt)  - * [(^)  « Xt(^)  t •••  j* 


so  that 


vS  H e»‘  .a ,2»t  + uVx*  / vS 


>2.2 


- + *“*  1 ^ « 


(q) 


^y  esqianding  u and  ths  algebraic  term  in  (o)  in  powers  of  s = 
and  combining  terms*  one  <*tains  finally  a relation 

-s^ 

e ® [ \/5  se®  erfo(s)  + 3s*’]  ds 


KA)  ~4 


r> 

J a 


e"®“  ds  - 


-.2 

X -O 


I 


J f ^ 2 

+ A / «‘®  15>^  + 2a*l  + 0(1/A*)  . (r) 
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The  error  in.vdL'»ad  in  extending  tte  integrals  of  (r)  to  s = «>  » ratteer 
than  8 * A f Is  of  order  e"*  ^ » and  so  does  not  change  tli©  asgmmtctlc 
expansion  indicated  iDy  (r}ii  evaluating  the  integrals  in  (r)»  one  finds 


1(A)  2 [1  - ^ 

A > 


(s) 


• 'f  ■ 


..  ^ f’  M:£Lsx— 


7 

+ • • • ^ 


t0:  sst&fsl^sAii 


-ISm,,  fffi-  Re(r),  Re(a}  > 0, 


Co) 


(d) 
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0(»slddr  a typical  integral  in  (c)®  for  e3n!i^:!le  that  appearing  in  the 
ooefficient  of  vithin  the  second  bracket t 

I = r . (b) 

L (1.177s 

If  ti^  e^onent  of  the  In  ^ factor  vere  instead  then  for 

He  (s)  > Ot  eq*  (d)  vould  give  fca*  I 

I (1“v)  (In  dv  ~ [1  ••  (^)  1 r(fl)9  (f^ 

'Jq 

But  it  is  readily  verified  that  bbth  sides  of  eq*  (f ) are  regular  funetlox^ 
of  the  ootoplex  variable  s for  BeCs)  > -I;  the  singularity  at  s = 0 
being  only  apparent.  Therefore » by  the  theory  of  analytie  continuation » 
the  equality  (f)  renaains  valid  for  ReCs)  > -I.  In  particular*  for  s = - 1/2 » 
(f)  gives 

I = [1  - 2^^]  r(-l/2),  (g) 

The  other  integrals  appearing  in  (c)  may  be  similarly  evaluted* 

Rram  eqs*  (o)»  (d)  one  thus  <^taina 


'Jo  \/r(w)  - u(v) 

-/t  V ^ r(^)  + r(|)  - 1 r(“  |)  (i  - 2^/^)  a^l 

+ ^ r(|)  - r(-  + I r(-  |)(i  - 2*2^/^+ 

- I r(-|)  (1  - 3’/^)  a^l  + 

utaioh  reducer,  to  eq.  17(46)  on  evaluating  the  gu&v^i  functions. 
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Aaroptotlo  Itese  Integral  IV(59). 

Elf  dlffd3rentiating  eqe  IV (58) 


4 X 


1 + 


T75 


+ *••  + 


b-  ln(tt*^^) 

^ 4 b.  ^ 

ir?75"'’6  «-o . 


u > u- , 


(a) 


&Tid  substituting  into  the  tanpsraturs  integs’al  r/(59)  there  results  after 
a elftange  of  variable  j 

^ C'  siisXj^  / — ,^4£. 


2_A 


i7S 


^ ge: 


1— ^ 


\27o 


“-iJg 

V 


1 


3 ..  V 


^ I 3 {„^j575 


'7 


-.1  [in  V + ln(u-u^)  - 2l 
» o 


U--U^ 

L.  O 


(b) 


Tna  tenss  In  (b)«  except  fof  the  last  atie,  yield  hy^rgeoinetric*' 
type  integrals  (inccssplete  Beta  functions) 

1 s ^ 

I (x)  = f ~ A [^(i.  s;  s+1;  1)  * X®  F(J?  S|  s+l*  x)]>  (o) 

® -'x  [/n-r  ® ^ 


r4)  r(e)  - 

« .,..1,— JL  5,8 

1)  s 


as 


X ■*  0^, 


(a) 


r(B  + j) 

valid  for  0<x<  1f  ReCs)  > 0.  Hence » by  the  theory  of  enalytle 
eontinuationii  (c)  is  valid  provided  only  that  0<x<  1»  8^0$  **'!»  *‘2;  ««<  » 
and  (d)  holds  ehen  ReCs)  '=■!•  0 < x < ci^  For  8 = 0,  (d)  is  meaningless, 
but  my  be  defined  by  a luting  procedure.  By  differentiating  (d)  vith 
respect  to  s at  s = - ^ , one  raadily  finds 


/. 


^ ^ " 2n  + [2  + In  x1  as  x *»  o"®"* 

x r'  v/W 


(e) 
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Thus  sqs*  ((l}t  (e)  gitra  for  (b) 


after  e rearrangsaisrits  ssvsrsl  tgrsss  «aj;oalllngs  ccaTOrls(m  with  eqs  (a)? 
the  last  grcup  of  terns  In  (f)  is  siiqply 


b^  - /<.z(u.  )1  = 


1 


since  by  assuuptlon  z(u^)  = 1»  Eq*  (f)  reduces  to  eq«  IV(?9)  of  the 
text  upon  evaluation  of  the  ganna  Amotions* 


0 
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♦ 


s 


4 


Ss  ^.laiarleal  Solution  tar  -bhe  Collapsinp  Bubble. 
The  system  IV (78)  may  be  uritteu 


(a) 

(b) 
(o) 
(d) 

(^) 


In  order  to  obtain  a sehemo  for  ntsaerioal  integrfl.tion,  subdivide  the  range 
of  values  of  u into  icitervals  defined  by  the  points 

0 = < •••  < Uj^  < < *•»  . it) 


The  intervals  eerrespoMing  to  (f)  are  in  general  not  equal;  but  chosen 
as  the  integration  proceeds.  For  oonvenienoe  wlte 

vi  ■ ^ ® <*' 


and  aasuma  that  f k ” ^ ^ 

are  known  for  0 < k < n. 

If  the  intei  val  - Uj^  is  aofflcientiy  shoxt;  0*  (u)  wltiiirt 

the  interval  may  bo  approximated  hy 


k vi  • " 


The  intee^’al  in  (c)  evaluated  at  the  point  u = moy  then  be  estimated  as 

r“n+i  ^ .ft 

/ iy)  'Ai-+i  - V dv  *=  > V,.  / 

'/ft  k=o 


= f ^ - ‘Vi  - Vi 


(1) 


y 


t 


Dsfisa 
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1^  = 0, 


^ - K+1  - 


3”5  B 


(J) 


o 


inan  accsrdlcg  to  eqs,  (i}y  (j)  and  (g}»  eq.  (c)  at  u " beocanss 

?n+1  + 3;^f  - ®n>' 


(k) 


Tho  valu©  of  ^ in  (k)  xaay  be  estiim-ed  In  tsrias  of  ty  as 

o:q)ansion  of  about  a value  © ~ § near  which  luieis  equllibriusx 

vapor  density  data,  say 


fn+1  = ? Vi)  = <5  [1  + 3,(e„^,  - e)  + 5^(8^,  - %)^  + (i) 

Thus,  for  the  first  few  steps  of  integration,  6-0,  suid  so  ^ ~ ^ 
Initially,  The  temperature  integral  relati^  (k)  becoaies 


= 1 


+ I + 
n 


[ ••  S)  - (6„  ~ 6)1  f 


(ni) 


la  which  tho  oaly  unknowns  are  and  easily 

solved  fui'  by  «n  Itei’etioii  prucedui'e  based  ou  an  alternative  form 

of  (m), 


V)  - ®)  ■= 


■ , Vlf 

^ - Vi  f 15,  + aj(e„^,  - Si  + - 1 


(n) 


To  obtain  , the  uifforonLiJil  equation  must  ba  usod»  At  oaoh 

*»•  f 

point  u « moke  tho  approodmations 
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W+l  - "n  J 

' V«h-vi/ 

2 <Vi  -Vi> 


sUu  ) K t 


“'A- 


= *A  » 


(h£Ll^]  + /'n' Vl\ 

,V  Vi  ~\i  W ' Vi  i 


> (o) 


„ J,  mttfu. 
O “«.VV 


^o««  «li««  \ 

2 “n'“n+r‘“n-n-r 


With  these  approacimatlons»  the  fornula 


z(ti)  = z +(ia“u)a*  + 
' ' n ' n n 


<"-V 


s« 

n 


is  exact  for  u = * Whoa  (u^  -•  ) equals  h = 

eqs.  (o)  give 


Vi  - V’ 


z” 

n 


,,  =!B±LliiL 

n 2h 


& .,  - 2b^  + z„  , 

- ,ati Is=I 


(p) 


BO  that  the  differential  equation  (a)  at  u = 1®  appraxlniatod  by  the 

dlffdronoe  equation 

n+1  ^7  n n~1 ' n+1 

* r < -1  • “ vi ' " ^ * **=>]  = “• 

n n 

Cdven  eq.  (q)  may  be  solved  for  the  positive  root 
being  known  i^om  equillbriun  vapor  pressuxt)  data. 
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t 


In  to  kssp  the  diffsrsaess  la  s;  s:;d  © sasall  aad  euaUre 

that  a poBltiTa  root  of  the  difference  equatim  (q)  sidata#  it 

beocioBa  nacasssuqr  to  deoreaee  h as  the  nuaarieal  integraticia  proeeeda* 
At  the  step  uhara  a new  value  of  h la  intr^ueedf  the  approcidiBatians 
(o)  rather  than  (p)  nust  be  '.ised»  The  differenee  equation  which  applies 
at  that  step  is  therefore  not  (q)»  but  ooe  obtainable  trm  (o)* 

8ta3rt  the  integration  * a fietitloiis  point  = - tij 
is  used*  Corresponding  to  the  initial  cenditicn  z’(0)  ss  0 end  the 
approadmation  (p)  $ one  Ksist  then  choose  z ^ = e.  • Since  a “ i ^ 

"•II  o 

difference  e<pmtion  (q)  for  n = 0 eiapHMes  to  the  linear  aquation 

«,  = 1 - ^ (1  + d,)  W,  s ^(e„)  s 0(D)).  (r) 


The  tenqperature  equation  loesses 


«i  = 


(•) 


z^(d^  '+d2  6^+»**) 


for  a = 0,  since  = § = 0,  ™ * 1,  For  sufficiently  aaiall  h, 

oq»  (s)  my  bo  approximted  by 


6 


i 


!iJLL 


3’r  ®1 


(t) 


D 9 P 9 


depend  on  0^  and 


It  tihould  be  noteil  that  t , d-#  o^i 

J I *• 

change  with  eftoh  new  oxiAnslon  of  ^ (ft)«  Caoause  these  pttraffletoi^5  j 

a«  t«H  os  b»  my  bo  oonetant  over  several  iitopa  of  integration » we  fitve 
not  given  mmiu  iiuui.(,M(s  which  depend  cn  u (i«.e»  ofi  n). 
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